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Abstract 

We study physics potential of Very Long Base-Line (VLBL) Neutrino-Oscillation 
Experiments with the High Intensity Proton Accelerator (HIPA), which will be com¬ 
pleted by the year 2007 in Tokai-village, Japan, as a joint project of KEK and JAERI 
(Japan Atomic Energy Research Institute). The HIPA 50 GeV proton beam will de¬ 
liver neutrino beams of a few GeV range with the intensity about two orders of 
magnitude higher than the present KEK beam for K2K experiment. As a sequel 
to the proposed HIPA-to-Super-Kamiokande experiment, we study impacts of ex¬ 
periments with a 100 kton-level detector and the base-line length of a few-thousand 
km. The pulsed narrow-band beams (NBB) allow us to measure the 
transition probability and the survival probability through counting experiments 
at large water-Cerenkov detector. We study sensitivity of such experiments to the 
neutrino mass hierarchy, the mass-squared differences, the three angles, and one 
GP phase of the three-generation lepton-flavor-mixing matrix. We find that experi¬ 
ments at a distance between 1,000 and 2,000 km can determine the sign of the larger 
mass-squared difference (ml — mf) if the mixing between and (the heaviest- 
or-lightest neutrino) is not too small; 2\UeJf‘{l — |Ke 3 p) ^0.03. The CP phase can 
be constrained if the iLesI element is sufficiently large, 2 |t 43 p(l — \Ueff‘) ^0.06, 
and if the smaller mass-squared difference (ml — m\) and the Ue 2 element are in the 
prefered range of the large-mixing-angle solution of the solar-neutrino deficit. The 
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magunitude \m\ — m\\ and the matrix element 17^3 can be precisely measured, but 
we find little sensitivity to — m\ and the matrix element C/e 2 - 
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1 Introduction 


Many neutrino experiments [0]-[lO| strongly suggest that there are flavor mixings in the 


lepton sector, and that neutrinos are massive. According to the atmospheric-neutrino ob¬ 
servations [^, the lepton-flavor-mixing matrix (the Maki-Nakagawa-Sakata (MNS) matrix 
pn| ) has a large mixing angle. Especially, the Super-Kamiokande (SK) collaboration 
reported that oscillates into the other species with maximal mixing. The K2K ex¬ 
periment, the current long-base-line (LBL) neutrino oscillation experiment from KEK to 
SK with L = 250 km and {Ey) = 1.3 GeV, obtained results which are consistent with 
the neutrino oscillation in the atmospheric-neutrino anomaly with 5m? ~ 3 x 10“^ and 
sin^ 26 r'U 1. The two reactor neutrino experiments CHOOZ and Palo Verde re¬ 
ported that no oscillation is found from and they exclude significant mixings between 
Pg and the other neutrinos. An important conclusion from these observations is that the 
atmospheric neutrino oscillation cannot be due to oscillation. Recently, the SK col¬ 
laboration showed an evidence that v,, oscillates into an active neutrino rather than sterile 
neutrinos 0. According to the results of solar-neutrino observations [0, the z/g flux from 
the sun is less than that of the prediction of the standard solar model 0 and the reduc¬ 
tion factor depends on neutrino energies. The most convincing explanation for this deficit 
is the z/g oscillation to the other neutrinos. Four possible scenarios of the solar-neutrino 
oscillation have been identified : the MSW BUI large-mixing-angle (LMA) solution, 
the MSW small-mixing-angle (SMA) solution, the vacuum oscillation (VO) solution 
and the MSW \ow-5m? (LOW) solution. Recently, the SK collaboration reported the en¬ 
ergy spectrum and the day-night asymmetry data and showed that the LMA solution is 
more favorable than the other scenarios |^. The SNO experiment which observes 


the solar neutrino flux with heavy water, showed conclusively, when combined with the 
SK flux data, that z/g oscillates into the other active neutrinos m- A consistent picture 
of three-neutrino oscillations with two large mixing angles and two hierarchically different 
mass-squared differences emerges from those observations, with the exception of the LSND 
experiment |]T^ which may indicate the existence of the fourth and non-standard (sterile) 
neutrino. 


Several LBL neutrino oscillation experiments 16, 17, 18, 19| and a short-base-line exper¬ 


iment [201 have been proposed to confirm the results of these experiments and to measure 


the neutrino oscillation parameters more precisely. The MINOS experiments |]^, from Fer- 
milab to the Soudan mine, with the base-line length of L = 730 km and (Ey) = 3.5 GeV, 
will start producing data in 2005. The observation of the survival probability Py^^y^ will 
allow us to measure the larger mass-squared difference and the mixing angle with about 
10% accuracy. Two LBL experiments, IGARUS and OPERA |jT^, from GERN to 
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Gran-Sasso with the base-line length of L = 730 km and at higher energies with {E^) ~ 20 
GeV have been proposed, and they may begin operation in 2005. The GERN experi¬ 
ments expect to observe the Vr appearance. Physics discover potential of these LBL 


experiments have been studied extensively K 


In Japan, as a sequel to the K2K experiment, a new LBL neutrino oscillation experiment 
between the High Intensity Proton Accelerator (HIPA) and SK has been proposed |]^ . 
The facility, HIPA, has a 50 GeV proton accelerator, which will be completed by the year 
2007 in the site of JAERI (Japan Atomic Energy Research Institute), as the joint project 
of KEK and JAERI. The proton beam of HIPA can deliver high intensity neutrino beams 
in the ~ 1 GeV range, whose intensity is two-orders-of-magnitude higher than that of 
the KEK 12 GeV proton synchrotron beam for the K2K experiment. The HIPA-to-SK 
experiment with L = 295 km and {E^) = 1.3 GeV will measure the larger mass-squared 
difference with 3 % accuracy and the mixing angle at about 1 % accuracy. 

All these experiments use conventional neutrino beams, which are made from decaying 
pions and Kaons that are produced by high-energy proton beams. The possibility of a 
neutrino factory has been discussed as the next generation of LBL neutrino-oscillation 
experiments 0. Here the neutrino beam is delivered from decaying muons in a muon- 
storage ring, where the stored muon energy may be in the several 10 GeV range. A neutrino 
factory can deliver very high intensity neutrino beams that are consist of the same amount 
of and (z7^ and Ve from /i+) with precisely known spectra. Possibility of constructing 
a neutrino factory in the JAERI site by upgrading the HIPA is now being extensively 
studied . 

In this paper we examine an alternative possibility of using conventional neutrino beams 
from the HIPA for Very-Long-Base-Line (VLBL) neutrino oscillation experiments, whose 
base-line length exceeding a thousand km A possible 100 kton detector in Beijing 


27| can be a target at about L = 2 ,100 km away. Physics capability of such experiments 


should be seriously studied because a neutrino factory may turn out to be too difficult or 
too expensive to realize in the near future. In order to take full advantage of conventional 
u beams, we examine the case of using pulsed narrow-band beams (NBB) and as a 
target we consider a large water-Cerenkov detector a la SK which is capable of measuring 
the z/^-to-z/g transition probability and the survival probability. We study sensitivity of 
such experiments to the signs and the magnitudes of the two mass-squared differences, the 
three angles and one GP phase of the three-flavor MNS matrix. 

This article is organized as follows. In section 2, we fix our notation and review the 
present status of the neutrino-oscillation experiments in the three-neutrino model. In sec¬ 
tion 3, we study the properties of the narrow-band neutrino beams that can be delivered 
by the HIPA 50 GeV proton synchrotron. In section 4, we study the signals, the back- 
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grounds and systematic errors of the VLBL experiments and present our findings for the 
prospects of experiments at the base-line length of 2,100 km and 1,200 km. Our results 
are summarized in section 5. 

2 Neutrino oscillation in the three-neutrino model 

In this section, we give the dehnition and useful parameterization of the 3x3 Maki- 
Nakagawa-Sakata (MNS) lepton-flavor-mixing matrix [^, and give constraints on its ma¬ 
trix elements and the neutrino mass-squared differences. 

2.1 The MNS matrix 

The MNS matrix is dehned analogously to the CKM matrix through the charged- 
current (CC) weak interactions, where the charged-current can be expressed as 

Jcc = 7^(1 -75)(M,c,t)^ + (e,7Z,r)7'^(l 1^2,(2.1) 

Here {i = 1,2,3) denotes the neutrino mass-eigenstates. The flavor-eigenstates of the 
neutrinos are then expressed as 


3 

2 = 1 


( 2 . 2 ) 


where a = e, fi,T are the lepton-flavor indices. 

The 3x3 MNS matrix has three mixing angles 
the following parameterization 


and three phases in general. We adopt 
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where V is the diagonal phase matrix with two Majorana phases, ip 2 and ip^. The matrix 
U, which has three mixing angles and one phase, can be parameterized in the same way as 
the CKM matrix. Because the present neutrino oscillation experiments constrain directly 
the elements, Ue 2 , Ues, and we hnd it most convenient to adopt the parameterization 
ii where these three matrix elements in the upper-right corner of the U matrix are the 
independent parameters. Without losing generality, we can take Ue 2 and to be real 
and non-negative. By allowing f/gs to have the complex phase 


Ue2, U ^3 > 0, Ue3 = [f/esl 6 


(0 < ^mns < 27r), 


(2.4) 


3 








the four independent parameters are Ue 2 , \ Ue 3 \ and ^^ns- other matrix elements 

of the U are then determined by the unitary conditions : 


Uel = 

U,l = 

Url = 


7l-|f/,3|2-|f/e2P, 

_Ue2Ur3±Uj^3U^ 

l-\Ue3\^ 

Ue2U^3 - Ur3UelU:^ 

l-lt/esP 


Ur3 = 

U,2 = 

Ur2 = 


UelUr3 - U,sUe2U:, 

l-lt/esP 

U^3Uel + Ue2Ur3Ut^ 

l-\Ue3? 


(2.5a) 


(2.5b) 


(2.5c) 


In this phase convention, f/gi, Ue 2 , and U-r 3 are all real and non-negative numbers, and 
the other hve elements are complex numbers. 

The Jarlskog parameter ||^ of the MNS matrix is defined as 


'^MNS 


Im (nv;vffitc,) = Im {u^,u;jjf,u:^) = (t/.,) , ( 2 , 6 ) 


where (a, {3) = (e, /r), (/i, r), (r, e) and (i, j) = (1, 2), (2, 3), (3,1). The last expression above 
is obtained in our phase convention. The two Majorana phases (^2 and do not contribute 
to the Jarlskog parameter. 

2.2 Constraints on the MNS matrix and the mass-squared dif¬ 
ferences 

The probability of finding the flavor-eigenstate (3 from the original flavor-eigenstate a at 
the base-line length L in the vacuum is given by 


F, 




m. 


(F^nr) 3 i 


i=i 


MNS/JO 




—7-7-* 

a3 


= - 4Re siti^ ^ sin" ^ sin" ^ 


2 Im 


UaiU^iUp 2 U^ [sin Ai 2 sin A 23 + sin A 31 ], 


( 2 . 


where Aij is 


Aij = -r^L ~ 2.534 


2 F. 


Jm|(eV 2 ) 

F,(GeV) 


■L(km), 


( 2 . 8 ) 
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with the neutrino energy Ey and the mass-squared differences Srnfj = — rrif. In partic¬ 

ular, the survival probability in the vacuum is 


= 1 - 4 {|t/„if/«2l^sin2 ^ ^ . (2.9) 

The oscillation probabilities of anti-neutrinos in the vacuum are obtained from those of 
neutrinos simply by reversing the sign of the Jarlskog parameter; 


P- - = P 

Ua^Va ) 




( 2 . 10 a) 

( 2 . 10 b) 


For instance, 


Py 


= -4Re 


I sin^ ^ + U,^U:,UaUl, sin 


2 ^23 


+ T 2 J^j^g {sin Ai 2 -|-sin A 23 -|-sin A 31 } . ( 2 . 11 ) 


The unitarity leads to the identities 

y~! ^ Pug^Ua = ^ P^a^Vg = ^ Pug^Ua ~ ^ ’ (2-12) 

f3=e,g,T f3=e,g,T fS=e,g,T (}=e,g,T 

The following approximations for the oscillation probabilities are useful in our study. 
When |A^ 2 l ^ I^i 3 l ~ 1) fhe oscillation probabilities can be expressed as 


Py 




2 ^13 


-|- |2Re (^^ 12 ^/ 32 ^/ 33 ^ 03 ) sin A ^3 ± 4 J]yj ^3 sin 


2 ^13 


A 


12 


+Re cos a, 3 ) A ?2 + 0(AL). (2.13) 

When 1 IA 12 I ~ 1 ^ 13 1; we may take into account finite resolution of A^j (A^^ 

A-j ± SA-j), and find 


Py 


^ A.,+iA.. 


2(5A 


27 A . . —5A . . 

J IJ ZJ 


dA^jPy^^y^ (A-^.) 
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,0T? fri TT* TT TT* A sin 5Ai2 siii 

+2Re < upiUP 2 Ua 2 cos ^12 + ^ 02 ^/ 32 ^/ 33 ^a 3 COS A 23 

, rr rr* rr tt* a sin^Agil 

+ ^a3^/33^/31^Ql cos A 33 j 


sin(5Ai2 , . ^ sin(5A23 , . ^ sin(5A3i] 

-hsmA^s^TT-hsmA33—-- ). (2.14) 

0 A 31 J 


sinozA^^ 

±2J„„, |s,„ + sm A,3^ 

In LBL experiments the uncertainty in dictates A^^, 


.A A 

rlA -. ~ —A- - 
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and the following two cases are relevant. When 

|iA,3| ~ 


A 
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E„ 


< 1 < 
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13 


E,. 


~ I^A^gl ~ I^A 


231 ) 


eq. (| 2 . 14 j) can be expressed as 

^ - 2 Re {u^,u;,u^2u:2 + u^2UhUp,u:, + u 

+ 2 Re \u^lU1^2^a2 cos A;^2} i 2 J^j^g {sin A32} 
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On the other hand, when 
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eq. (|2.14j) is simplihed to 
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Note that in eqs. (|2.13|) , (|2.14|) , and (|2.17|) , stands for Im PQiP^iP^ 2 ^a 2 

eq. 


; see 
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All the above formulas remain valid for the neutrino oscillation probabilities in the 
matter, by replacing the mass-squared differences and the MNS matrix elements with the 
effective ones in the matter, 


A 




^ij 7 Ufy 




J., 


J., 


( 2 . 20 ) 


as long as the matter density remains the same along the base-line. The dehnitions of 
the effective parameters Ajj and Uai are given in section and is obtained from 
eq.(|2.6|) by replacing all Uai^s by Uai^s. 

In the following, we summarize the constraints on the neutrino mass-squared differ¬ 
ences and the MNS matrix elements from the recent neutrino-oscillation experiments; the 
atmospheric-neutrino anomaly |^, the CHOOZ reactor experiment |^, and the solar- 
neutrino deficit observations []^, jTO 


2.2.1 Atmospheric-neutrino anomaly 

A recent analysis of the atmospheric-neutrino data from the Super-Kamiokande (SK) ex¬ 
periment 10 finds 

0.88 < sin2 20^^^ <1.0, (2.21a) 

1.6 X 10"=^ < (eV^) < 4.0 X 10“^ , (2.21b) 

from the —>■ survival probability in the two-flavor oscillation model; 


P. 


= 1 


sm 


20ATM sin 


' Sm? ^ 

ATM 

4K. . 


( 2 . 22 ) 


The base-line of this observation is less than about 10^ km for the Earth diameter, and 


the typical neutrino energy is one to a few GeV. The survival probability eq. (p.9|) may then 
be expanded as 


=1-4 


U. 


/i3 


1 - 


u, 


/x3 


. 2 ^13 
sm - 


+ 2 


U, 


/i2 


u. 


fi3 


AusmAi3 + 0{Al^). (2.23) 


When IA 12 I -C 1, we may neglect terms of order A 12 and obtain the following identihcation 


4\U,s\\l-\U,3\^)=sm^29, 


5ml^ 


= 5rrr. 


(2.24a) 

(2.24b) 
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The independent parameter (> 0 in our convention) is then 




sin^ 29 


ATM 


/ 


72 . 


(2.25) 


The magnitude of the neglected terms in the above approximation is largest when the 
large-mixing-angle solution of the solar-neutrino dehcit is taken. For = 10“^ eV^, 

Ey = 1 GeV, and L = 10^ km, we have IA 12 I — SrnlSL/2E,) 1, and a more careful 
analyses in the three-neutrino model are required to constrain the model parameters. The 
results of such analyses [^] show that the identihcations eqs. (p.24|) remain valid approxi¬ 
mately even for the large-mixing-angle solution. 


2.2.2 Reactor neutrino experiments 

The CHOOZ experiment [Q measured the survival probability of z7g, 


P- - = 


sm 


^^CHOOZ 




(2.26) 


and it was found that 



' 0.10 

for 

dm? 

CHOOZ 

> 3.5 X 10-3 eV^ 

sin 20 (^jjqq 2 ^ 

0.18 

for 

CHOOZ 

> 2.0 X 10-3 eV^ 


0.52 

for 

CHOOZ 

> 1.0 X 10-3 eV^ 


(2.27) 


The base-line length of this experiment is about 1 km and the typical anti-neutrino en¬ 
ergy is 1 MeV. For those L and Ey, the Earth matter effects are negligible, and IA 12 I = 
5m?^oi^{L/2Ey) can be safely neglected even for the large-mixing-angle solution. The sur¬ 
vival probability of the three-neutrino model is then approximated by 


Pv.^v. = 1-4 \U„\^ (1 - \Uj) si 
and we obtain the identihcations : 


sm 


Ai3 


+ 0(A 


12 


\5m\.^\ = 5m 


CHOOZ ’ 
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CHOOZ* 


(2.28) 


(2.29a) 

(2.29b) 


With the above identihcations, we hnd that the element jf/ggl eq. (|2.29a|) is constrained 
by eq.( 2.27 ) in the region of allowed by the atmospheric-neutrino oscillation data 














through eq. (|2.21b|) and eq. (p.24b|) . The independent parameter is now constrained by 
eq. (|2.27|) through the identihcation 





- sin^ 29 


CHOOZ 


/ 


72. 


(2.30) 


2.2.3 Solar-neutrino deficit 


Dehcit of the solar neutrinos observed at several terrestrial experiments il, [10| have been 
successfully interpreted in terms of the z/g — i>x {^x 7^ ^e) oscillation 


f5rn? \ 

= 1 - sin^ 293„ sin^ ( ^ 

in the following four scenarios 0, P, [T^ . 

MSW large-mixing-angle solution (LMA) : 

0.7 < sin2 20goj^ <0.9, (2.32a) 

3 X 10"® < < 15 X 10"®. (2.32b) 


MSW small-mixing-angle solution (SMA) : 

1.2 X 10"=^ < sin2 20goj^ <12x10-7 (2.33a) 

0.3 X lO"*^ < 5m77eV7 < 1 X 10-®. (2.33b) 

MSW low-5m^ solution (LOW) : 

0.8 < sin^ 29 ^^^ < 1, (2.34a) 

10“® < <5m77eV7 < 2.5 X 10"^. (2.34b) 

Vacuum Oscillation solution (VO) : 

sin^ 29 ^^^ ~ 0.9 , (2.35a) 

5m77eV7~10-7 (2.35b) 


The SK collaboration reported that their data on the energy spectrum and the day-night 
asymmetry disfavor the SMA solution at 95% C.L.. Recently the SNO collaboration gives 
us the hrst direct indication of a non-electron and active flavor component in the solar 
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neutrino flux. Because S> 5m? the z/„ survival probability in the three-neutrino 

ATM SOL e t' J 

model can be expressed as 


Pu 


1 - 2 It/esi' (l - It/esi') - 4 \Uel\^ \Ue 2 f siu^ + O 


E„ 




(2.36) 


where terms of order [A 13 ^ in eq. (|2.14|) are safely neglected. The energy- 

independent dehcit factor, 2 |f/g 3 p(l — |17e3p); should be smaller than 5% by the CHOOZ 


constraint eq. (p.27|) if > 3.5 x 10“^eV^. Because we need only rough estimates of 

the allowed ranges of the MNS matrix elements^, we ignore the small energy-independent 
dehcit factor and interpret the results of the two-havor analysis eq. ( p.32a|) -eq. (p.35a)) by 
using the following identihcations 


SOL ) 




By using unitary condition, the independent parameter U ^2 is obtained as 


(2.37a) 

(2.37b) 


V2 = 


(1 - Vest - yoHfTiTTtaTT^) 


1/2 


(2.38) 


2.3 Neutrino mass hierarchy 

All the above constraints on the three-neutrino model parameters are obtained from the 
survival probabilities which are even-functions of 5m‘fy We have made the identihcation 


'^"^SOL = 


5m? 


12 


•C 


5m\^ 


= 5m~ 


(2.39) 


which is valid for all the four scenarios of the solar neutrino oscillation. 

There are four mass hierarchy cases corresponding to the sign of the 5mjj, as shown 
in Fig. 1^ and Table [^. We name them the neutrino mass hierarchy I, II, III and IV, 
respectively. 

If the MSW effect is relevant for the solar neutrino oscillation, then the neutrino mass 
hierarchy cases II and IV are not favored, especially for the LMA and SMA solutions. The 
hierarchy I may be called ‘normal’ and the hierarchy III may be called ‘inverted’. Within 
the three-neutrino model, there is an indication that the normal hierarchy I is favored 
against the inverted one III from the Super-Nova 1987A observation [^. Nevertheless a 
terrestrial experiment is needed to determine the neutrino mass hierarchy. 


^See for example, more detail discussions in 
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Table 1: The four neutrino mass hierarchies and the corresponding sign assignments for 

(5m^2 6ml^. 
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Figure 1: Schematical view of the four-types of neutrino mass hierarchy. 


We notice here that there are two types of mass eigenstate for the neutrinos. The states 
Vi with the mass rrii {i = 1, 2, 3) appear in the definition eq.( p.2|) of the MNS matrix, whose 
elements are constrained by the existing experiments that measure essentially the neutrino- 
flavor survival probabilities. Since the survival probabilities eq.( p.9| ) do not depend on the 
sign of the mass-squared differences, these constraints do not depend on the neutrino mass 
hierarchy. Because the MNS matrix elements are constrained uniquely by the neutrino- 
flavor survival probabilities, we may call these states as ‘current-based’ mass-eigenstates. 
We find this basis most convenient for our study in this paper. On the other hand, the 
‘mass-ordered’ mass-eigenstates v[, whose masses satisfy 


rrii < m2 < rUg , 


(2.40) 


are useful when studying the high-energy behavior of the neutrino mass matrix |Q, the 
matter effects on the neutrino-flavor oscillation [T^, |T^, and when studying the lepton- 
number violation effects which are proportional to the magnitudes of the Majorana masses. 
The relation between the current-eigenstates and the two mass-eigenstates is 




( Vi ^ 
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(2.41) 


where V is the MNS matrix in the mass-ordered mass-eigenstate base. It can be obtained 
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relate the two mass-eigenstates 



= 


^2 
v'-i 
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for the neutrino mass hierarchy I, II, III and IV, respectively. 

2.4 Neutrino oscillation in the Earth matter 

Neutrino-flavor oscillation inside matter is governed by the Schrodinger equation 


d I 


dt 


V,, = 


V 


2E,, 


Ho + 


a 0 0 
0 0 0 
V 0 0 0 


/ 


J V 


/ 


= H 


where Hq is the Hamiltonian in the vacuum, 

/ 0 0 

Ho = U 0 6m 

V 0 0 

and a is the matter effect term jl^. 


0 

0 

5ml^ 




= 2V2GFneE^ = 7.56 x 10-^(eV^ 


P 


^g/cm'^ 


E^ \ 
GevJ ■ 


(2.42) 


(2.43) 


(2.44) 


(2.45) 


(2.46) 


(2.47) 


Here rig is the electron density of the matter, Gf is the Fermi constant, and p is the matter 
density. In our analysis, we assume that the density of the earth’s crust relevant for the 
VLBL experiments up to about 2,000 km is a constant^, p = 3, with an uncertainty of 
Ap = 0.1; 


p (g/cm^) = 3.0 ±0.1 


(2.48) 


^more detail discussion of the matter profile in Ref. [|35| . 
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The Hamiltonian in the matter H is diagonalized as 


H 


. / Ai 0 

0 As 
\ 0 0 


0 

0 

As 


U\ 


(2.49) 


by the MNS matrix in the matter U. The neutrino-flavor oscillation probabilities in the 
matter 






(2.50) 


takes the same form as those in the vacuum eg. ( |2.7|) , where the elements are replaced 
by and the terms A-j are replaced by 



2Ei, 2Ey 


(2.51) 


In Fig. g we show the Pj^-dependence of the effective mass-squared differences and the 
MNS matrix elements inside of the matter at p = 3.0 g/cm^. The curves are obtained 
for \5m\.^\ = 3.5 x 10“^ eV^, \ 5 m\ 2 \ = 15 x 10“® eV^, sin^ 20^.^,^ = 1.0, sin^26'gQj^ = 0.8, 
sin^ 26 '(^jjqq 2 = O-li = 0°. For the mass hierarchy III and IV, < 0, and for 

II and IV, 5m\2 < 0. The elements affect the survival probability in the 

matter, whereas the terms Ufj.iU*^ affect the transition probability ; see eq. (|2.9|) and 

eq. o. whose expressions remain valid in the matter by replacing Uai and Aij by Uai and 
Aij, respectively. 

Between the hierarchy I and II (III and IV), the sign of the smaller mass squared 
difference 5m\2 is different. Even though the individual terms behave quite differently in 
the matter, we And that neither Py^^y^ nor Py^^y^ depend strongly on the sign of 
in the range of (1 ~ 6 GeV) and L (300 ~ 2,100 km) that we study. On the other 
hand, we And strong dependence of the transition probability Py^^y^ on the sign of 
between the hierarchy I and III (II and IV). We show in Fig. ^ the > z/g amplitude 

(2.52) 

in the complex plane for the same parameter set at four points. The three terms in 
the above sum are shown in the complex plane as two-vectors, whose sum is chosen to 
lie along the horizontal axis. The absolute value squared of the sum gives the probability 
Py^^y^. The left figures are for the normal hierarchy I and the right figures for the inverted 
hierarchy III. The double circle shows the origin, and the solid-circle (with solid line). 
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solid-square (with long-dased line), open-circle (with short-dashed line), open-square (with 
dotted line) shows the amplitude for = 0°, 90°, 180° and 270°, respectively. The four 
hgures in each column show the amplitudes at the same LjE^, = 420km/GeV; from the 
top to the bottom, we show the amplitudes in vacuum, in the earth crust (p = 3.0 g/cm^) 
at L = 295 km, L = 1, 200 km, and L = 2 ,100 km. We can clearly see from these hgures 
that the to transition amplitude increases by the matter effect at higher energies 
and hence at a larger distance in case of the normal hierarchy I, whereas the amplitude 
decreases signihcantly in case of the inverted hierarchy III. The difference in the transition 
probability is more striking after the amplitude is squared. 





Figure 2: Mass-squared differences and the MNS matrix elements inside of the matter 
(p =3.0 g/cm^) for the four neutrino mass hierarchies, I (solid lines), II (long-dashed lines), 
III (short-dashed lines) and IV (dotted lines). The curves are obtained for \Sml^\ = 3.5 x 
10 “^ (eV^), |5m^2l = 15 X 10 “^ (eV^), sin^ 26'^.j,j^ = 1 . 0 , sin^ 26 ^ 3 ( 3 ^ = 0 . 8 , and siv? 29 = 

0 . 1 . 

Before closing this section, we give a useful relationship between the z/q, —z/p transition 
and the z/q, —z7p transition which may be valid in the terrestrial LBL experiments where 
both the accelerator and the detectors are near the earth surface. The oscillation of anti¬ 
neutrinos in matter is given by the Schrodinger eqnation: 




/ —a 

0 

Ml 

Ho + 

0 

0 

0 


1 0 

0 

oyj 




^ W \ 


= H 


\jyr ) 


\ ) 


(2.53) 


The Hamiltonian in the vacnum is identical to the one governing the neutrino-flavor oscil- 
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Figure 3: The v^. amplitude in the mater (p = 3.0 

g/cm^) at L/Ei, = 2,100 km/5 GeV for \Sml^\ = 3.5 x 10“^ eV^, |5mf2| = 15 x 10“® eV^, 
sin^ = 1.0, sin^ 26*gQL = 0.8, and siv? = 0.1. The three terms in the sum are 

shown in the complex plane, where the real part of the sum gives The four figures 

in the left are for the mass hierarchy I and the hgures in the right are for the hierarchy 
III. The solid-circle, solid-square, open-circle, open-square are for = 0°, 90°, 180° and 
270°, respectively. The double circle shows the origin. 
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lation eg . (12.461) , 


^o = ^£/t 


/ 0 0 

0 5m2i 

V 0 0 


0 

0 

5ml^ 


U = Hi = H 


0 • 


(2.54) 


By comparing the total Hamiltonians H and H; 


rim 2 

Tj 


2aE 


+ \U. 


e3 


LV 


u,3u:, 

Ur3U:, 


Uesu;, 

\U,3? 

Ur3U;, 


Uesu:, 

u,3u:, 

l^rSp 
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6ml 


13 


/ 


Ue2U;, 

Ue2U:2 \ 



u,2u:. 

\U,2\^ 

u,2u:. 


V 

Ur2U:2 

Ur2U;, 

\Ur3? ) 



, (2.55a) 


— ^6i7^ 
2E 



/ 

2aE 1 ITT |2 

+ l^esl 

u:sU,s 

u:sU,3 \ 



u;sUe3 

\U,3\^ 

u;,u^3 


1 

U:sUe3 

CO 

\Ur3? y 


dmls 


/ 

\Ue2? 

u:2u,2 

U:2Ur2 V 



U;2Ue2 

K2? 

U;2Ur2 

, (2.55b) 

V 

u:^u,2 

U*r,U,2 

\Ur3\^ )\ 



we find 




(2.56) 


Because the Hamiltonian —H governs the oscillation in the reversed time direction, we 
find that the following identities hold 




= Pj^ 


Ua — 


Ua —’ 




= 


Vql — 


Uoc —^ 

>I//3 

Pv"" 


= PJ^ 



>Uf3 



p IV 


= P ^ 


^ I'd— 

>VI3 


>I//3 


(2.57a) 

(2.57b) 

(2.57c) 

(2.57d) 


if the matter density along the baseline is symmetric under the reversal of the beam di¬ 
rection, i.e under the exchange of the injector and the detector. This condition is met 
approximately for all terrestrial LBL experiments where both the accelerator and the de¬ 
tector are on or near the earth surface. 

In the following, we therefore give results for all the four hierarchy patterns but only 
for the neutrino beam. Oscillation probabilities for anti-neutrino beams are then obtained 
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according to the rule eqs. (p.57|) , while the CC and neutral-current (NC) event rates are 
obtained after multiplying the ratio of the anti-neutrino and neutrino cross sections on the 
target. 


3 Narrow-Band Neutrino Beams with HIPA 

3.1 The KEK-JAERI joint project 

The KEK-JAERI joint project on HIPA is a proton accelerator complex and associ¬ 
ated experimental facilities which will be constructed in the site of JAERI, Tokai-village, 
60 km north-east of KEK. The project consists of 400 MeV Linac, 3 GeV and 50 GeV syn¬ 
chrotrons. The design parameters of 50 GeV machine are listed in Table ^ together with 
some other proton machines for LBL experiments. The intensity is 3.3x10^^ protons/pulse 
(PPP) cind the repetition rate is 0.275 Hz. The power reaches 0.75 MW which is 2 orders 
of magnitudes higher than the KEK 12 GeV proton synchrotron (PS). The accelerators 
in the facility will be in the power frontier in the world. The facility is approved by the 
Japanese government in December, 2000 and the construction will take 6 years from 2001. 

In the following discussion, 10^^ protons on target (POT) is adopted as a typical 1 year 
operation. This corresponds to about 100 days of operation with the design intensity. 

Table 2: Gomparison of accelerator parameters used for LBL projects 



Energy 

(GeV) 

Intensity 

(10l2ppp) 

Rep. rate 
(Hz) 

Power 

(MW) 

HIPA 

50 

330 

0.275 

0.75 

NuMI 

120 

40 

0.53 

0.41 

KEK-PS 

12 

6 

0.45 

0.0052 


3.2 Neutrino Beams for LBL experiment between HIPA and SK 


As a first stage neutrino experiment at this new facility, long baseline experiment from 
HIPA to SK has been planned and discussed seriously by the JHF neutrino working group 
m- Before going into the description of higher energy beam for VLBL experiment, we 
briefly introduce the beam for the LBL experiment. 

Major purposes of the HIPA-to-SK experiment are 1) precise measurement of oscilla¬ 
tion parameters in disappearance, 2) to discover z/g appearance. The principles of the 
experiment are 
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• Use of low-energy narrow band beam (NBB) whose peak energy is tuned at the 

oscillation maximum. Since the distance between HIPA and SK is 295 km, the peak 
energy should be around Ei, = 0.7 ~ 1.2 GeV for the region of allowed by the 

SK observation p. 

• Neutrino-energy is kinematically reconstructed event-by-event from the measured 
lepton momentum by assuming the charged-current quasi-elastic (CCqe) scattering. 
Inelastic scattering with invisible secondary hadrons mimic the CCqe interactions 
and smears the measurement. Below Ey 1 GeV, interaction is dominated by 
the CCqe interaction. So the low-energy beam with small high-energy tail is favorable 
for this method. 


Currently, three beam options are being considered, namely the wide-band beam (WBB), 
the NBB and the off-axis beam (CAB); 

WBB: Secondary charged pions from production target are focused by two electromagnetic 
horns [^. Since the momentum and angular acceptances of the horns are wide, 
resulting neutrino spectrum are also wide. The advantage of the WBB is the wide 
sensitivity in But backgrounds from inelastic scattering of neutrinos from high- 
energy neutrinos limit the precision of the oscillation-parameter measurements. 

NBB: Two electromagnetic horns have their axis displaced by about 10°, and an dipole 
magnet is placed between them to select pion momentum. Resulting spectrum has a 
sharp peak and much less high energy tail than WBB. 

CAB: The arrangement of beam optics is almost the same as WBB, i.e. coaxially aligned 
two horns. The axis of CAB is intensionally displaced from the SK direction by a few 
degrees. Pions with various momenta but with a hnite angle from the SK direction 
contribute to a narrow energy region in the spectrum . The CAB can produce a 
factor of 2 or 3 more intense beam than NBB. But unwanted high-energy component 
is larger than NBB. 


The length of the decay pipe is chosen to be relatively short, 80 m for all the conhgurations. 
This is because 1) high-energy neutrinos do not improve the measurements, 2) longer pipe 
costs very much due to the heavy shielding required by the Japanese radiation regulation. 


In [^, the WBB is used only in the early stage of the project in order to pin down 


at about ~ 10% accuracy^]. Typical expected spectra of those options are plotted in Fig. § 
and the flux and number of interactions are summarized in Table 

^ Very recently the JHF-SK neutrino working group has modified its strategy and the beam configura¬ 
tion. In the most recent plan, the OAB will be adapted for the LBL experiment and the decay pipe length 
will be 130 m. For more details see Ref.IM. 
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3.3 High Energy Narrow Band Beam for VLBL Experiments 

In order to explore the physics potential of the VLBL experiment with HIPA, we need 
to estimate the neutrino flux whose spectrum has a peak at higher energies. We study 
the profile of such beams at distances of 1,200 km and 2,100 km by using Monte Carlo 
simulation. In this subsection, we describe the beam in detail. 

First, we chose the decay pipe length to be 350 m. For the baseline length currently 
under consideration, 1,200 km to Seoul and 2,100 km to Beijing, oscillation maximum lies 
at F'y = 2 ~ 4 GeV and 4 ~ 7 GeV, respectively, for 6 ml^ = (2 ~ 4) x 10“^ eV^. In order 
to make 5 GeV neutrinos for example, we need pions of momentum about 10 GeV. The 
10 GeV pions run about 560 m during their life. Therefore we need a long decay pipe of 
several 100 m for efficient neutrino production. Gonsidering the site boundary of JAERI 
and layout of accelerators, maximum decay pipe length is about 350 m. 



Figure 4: Typical spectra of the GG interactions in the absence of neutrino oscillation 
in the HIPA-to-SK LBL experiment. The solid, dashed and dotted histograms are spectra 
for WBB, OAB(2°) and NBB {{Ptt) = 2 GeV), respectively. 


Secondly, we adopt the NBB conhguration. Use of WBB at these high energies has at 
least two disadvantages; 

(i). The reconstruction of neutrino energy is difficult. At multi GeV region, inter¬ 
actions are dominated by deep-inelastic scattering with multi-pion production. It 
is difficult for a water Cerenkov detector to make such measurements. It is a non¬ 
trivial exercise to construct a 100 kton-level detector of a reasonable cost which has 
the capability of reconstructing the neutrino energy at several GeV range. See e.g. a 
proposal in ref . 
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Beam 




^e-^tot 

NBB((p^) = 2 GeV) 

7.0 

870 

620 

6.8 

OAB (2°) 

19 

3100 

2200 

60 

WBB 

26 

7000 

5200 

78 


Table 3: Summary of beam simulation of the HIPA-to-SK LBL experiment. The neutrino 
flux $ at SK is in the unit of 10®/cm^/year, A"tot and Nqq are the number of total and CC 
interactions, respectively in SK’s hducial volume of 22.5kton for 1 year (lO^^POT) in the 
absence of neutrino oscillation. 


(ii). The construction of the beam line costs very much. In the WBB conhguration, the 
proton beam which passes through the production target goes through the decay 
pipe all the way down. The intensity is still of the order of 10^^ ppp even beyond 
the target. In order to shield the extremely high radiation, we need a considerable 
amount of shielding around the decay pipe. Constructing a decay pipe of several 100 
m with heavy shielding is unrealistic. 

The OAB also has the second disadvantage. Therefore we choose high energy NBB for the 
present study, since it does not suffer from the above disadvantages. With an ideal NBB, 
neutrino energy reconstruction is not necessarily be done by the detector. It is possible to 
design a NBB beam line where the 50 GeV proton beam does not enter the decay pipe. 
Simulation of high energy WBB is done only for comparison. 

For the purpose of the focusing secondary pions, we adopt the quadrupole (Q) magnets 
instead of horns. In general, focusing by Q magnets has smaller angular and momentum 
acceptance than the horn focusing. The reasons why we choose the Q focusing are 


• Q focusing gives narrower neutrino-energy spectrum 

• High energy pions of ~ 10 GeV are emitted at smaller angles and hence reasonable 
acceptance for those pions can be obtained by the Q optics. 

• The Q-magnet can be operated at low DG current of several kA. Gompared with the 
horn magnets which require pulsed operation with a few 100 kA, much more stable 
operation can be expected. 


We used GEANT 39 for the beam-line simulation to estimate the neutrino flux. A 


target, Q magnets, a dipole magnet, and a decay pipe are put into the geometry. The 
target is Gu rod of 1 cm diameter and 30 cm length. The length corresponds to about 
2 nuclear interaction length. GGALOR code is used for hadron production in the 
target. Every secondary particles produced in the target are tracked. The beam line optics 
assumed for the present study is drawn in Fig. The secondary pions from the target are 
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Dipole 

Magnet 


Q 



Target x 

100 cm 

z 

Figure 5: Beam line optics for high energy narrow band beam 

focused by the following 4 Q magnets and bent by 10° by a dipole magnet. The optics is 
not fully optimized. Several tens of % increase in flux could be expected by tuning the size 
or held of the Q and bending magnets, the position of the target etc. 

In Fig. P, some typical spectra obtained by the MC simulations are plotted. Spectra 
with a narrow peak structure are generated. The peak has a sharper edge in the high- 
energy side than in the low-energy side. The trailing edge in the low-energy side comes 
from pions with hnite angle. We observe a small secondary peak at an energy about twice 
the primary peak position. The second peak comes from Kaon decays. In the left-hand- 
side hgure, we show the hux of the neutrino at 2,100 km away from HIPA in units of 
10'^/400MeV/cm^/year, where 10^^ POT is assumed for one-year operation. Three types 
of the NBB where peak energy is at about 3, 5, 8 GeV, and the high-energy WBB spectra 
are shown. The right-side figures show the expected number of CC events par year, in 
the absence of neutrino oscillation, for a 100 kton detector. The cross section has been 
obtained by assuming that the target detector made of water. 

Spectra for each neutrino species in the 5 GeV and z7^ beams are shown in Fig. |^. 
The flux ratio of i^eO^e) to is 0.9% (0.8%) in total, and 0.2% (0.1%) at the peak 

energy for {y^) beam. In the left hgure, we show the z/^, z/g, z7^ and W spectrum of the 
NBB z/^ beam, while in the right hgure the corresponding spectrum are shown for the NBB 
beam. The number of wrong-sign z7^ (z/^) GG interactions in z/^ iy^) beam is about 
0.4% (4%) of the right-sign interactions. Although the hux of the z/^ beam and the 
beam are almost the same, as well as the fraction of contaminated neutrino hux, the z7^ 
beam suhers from 10 times higher wrong-sign events because of a factor of 3 smaller 
GG interactions than the z/^ GG interactions oh water target at these energies. 

The results of the simulations are summarized in Table ^ and Table In Table |^, we 
show the expected number of GG and NG events for the 3 GeV and 6 GeV NBB’s for 
100 kton-year (lO^^POT) at four typical distances, 300 km, 700km, 1,200 km, and 2,100km 
from HIPA. The upper numbers in each row and column show the numbers of events 
without oscillations, while the lower numbers are calculated by using the three neutrino 
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Figure 6: Typical spectra of narrow band beams. The left figure shows the flux and the 
right one gives the number of charged current interactions in the absence of oscillation 
at L = 2,100 km. The solid, dashed and dotted histograms are 3, 5 and 8 GeV NBB, 
respectively. The thin solid histograms show the WBB spectra for comparison. 




Figure 7: Composition of neutrino species for a typical NBB. Solid, dashed, dotted and 
dot-dashed histograms correspond to z/^, z/g, z7^ and z/g, respectively. The left and right 
hgures shows 5 GeV z/^ and z7^ beam, respectively. 
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Table 4: Expected number of interactions of NBB. Assumed parameters are hierarchy 
= 3.5 X 10-3 eV^, = 1.0, = 10 x lO"® eV^, sin^ 2030^^ = 0.8, 

sin^ 26 '(.jjqq2 = 0.06, = 0° and p = 3 g/cm^. First and second lines for each set of 

baseline length and the peak energy indicate the number of interactions without and with 
oscillation, respectively. 


L 

-^peak 

cc 

CC 

CC 

Fe CC 

Nnc 

300km 

3GeV 

7495. 

43.0 

55.0 

0.90 

2540.9 



5903. 

22.0 

105.0 

1.20 

2540.9 


6GeV 

13321. 

44.0 

82.0 

1.90 

4457.4 



12400. 

21.0 

110.0 

1.70 

4457.4 

700km 

3GeV 

1376. 

7.9 

10.1 

0.17 

466.7 



382. 

3.8 

43.9 

0.28 

466.7 


6GeV 

2446. 

8.1 

15.0 

0.35 

818.7 



1699. 

3.5 

40.6 

0.36 

818.7 

1,200km 

3GeV 

468. 

2.7 

3.4 

0.05 

158.8 



85. 

0.8 

18.1 

0.13 

158.8 


6GeV 

833. 

2.7 

5.1 

0.12 

278.6 



297. 

1.1 

25.3 

0.13 

278.6 

2,100km 

3GeV 

153. 

0.9 

1.1 

0.02 

51.9 



119. 

0.5 

2.4 

0.05 

51.9 


6GeV 

272. 

0.9 

1.6 

0.04 

91.0 



47. 

0.4 

13.2 

0.06 

91.0 


Table 5: Expected number of interactions of NBB. Assumed parameters are the same 
as in Table 


L 

'^peak 

CC 

17^ CC 

Ve CC 

Fe CC 

tV^NC 

300km 

3GeV 

160.0 

2871.0 

4.0 

12.7 

1184.7 



66.0 

2261.0 

4.9 

33.6 

1184.7 


6GeV 

141.0 

4076.0 

6.7 

13.3 

1584.8 



77.0 

3793.6 

5.4 

22.4 

1584.8 

700km 

3GeV 

29.4 

527.0 

0.7 

2.4 

217.6 



13.9 

136.0 

1.1 

16.6 

217.6 


6GeV 

25.9 

748.6 

1.2 

2.5 

291.1 



11.5 

518.5 

1.2 

10.8 

291.1 

1,200km 

3GeV 

0.0 

179.0 

0.2 

0.8 

74.0 



3.5 

28.0 

0.5 

7.0 

74.0 


6GeV 

00 

bo 

254.7 

0.4 

0.8 

99.1 



3.4 

88.1 

0.5 

7.5 

99.1 

2,100km 

3GeV 

3.2 

58.6 

0.1 

0.3 

24.2 



1.9 

47.6 

0.2 

0.8 

24.2 


6GeV 

2.8 

83.1 

0.1 

0.3 

32.3 



1.5 

12.9 

0.2 

4.1 

32.3 
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model for the following parameters : 


( 5m?__ 




= ( 3.5 X 10■^1.0 X 10■^ 1.0, 0.8, 0.06, 0° ) (3.1) 


with the neutrino mass hierarchy I and for a constant matter density of p = 3 g/cm^. Be¬ 
cause all the three neutrinos have identical NC interactions, the two numbers are identical 
in A^nc column. All the upper numbers simply follow the 1/L^ rule of the flux at a distance 
L. 


In Table 5, we show the corresponding numbers for the T7^ NBB’s. The number of 
expected events are about a factor of 3 smaller than the corresponding ones in Table 4 
because of the smaller CC and NC interactions of off nucleus target. 


Details of all the NBB’s generated for this study are available from ||^ 


3.4 Parameterization of the high-energy NBB 

In the numerical studies of the next section, we make the following parameterization of the 
NBB with a single peak at = i^peak- 


/ P \ ft(®peak) 1 / CCfrp 

( 3 . 2 ) 


where is running from 0 to A'peak; /(-^'peak) and 6(i?peak) are parameterized as 

/(Bp..k) = 3.3£;J„p - 76.8Bp..k + 520. (3.3a) 

5(£p..t) = 20.3B-i" + 2.8, (3.3b) 


where T^pgak is measured as unit of GeV. M = 100 kton stands for the mass of the 
detector, = 6.017 x 10^^ is the Avogadro number, ^{E^) is the flux (in units of 
/GeV/cm^/lO^^POT) at L = 2,100 km. The and are, respectively, the and Ug 
GG cross sections per nucleon off water target [^, and their ratio is approximately given 
by 


r 1.0 - 0.056E-0-48 {0.7 ^E^) 

crJ^(E,)/af^(E,) ~ i 0.83 + 0.16Ep (0.3 ^0.7) (3.4) 

[ 0.879(i7i. - m^)/(0.3 - m^) (m^ ^E^ ^ 0.3) 


where = 0.11 GeV represents the muon mass. This parameterization allows us to 

study the effects of changing the peak energy of the NBB continuously. We show in 
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Figure 8: The NBB’s from HIPA (histograms), and our parameterization (solid lines). The 
horizontal axis is the neutrino energy {E^^) and the vertical axis is the beam flux times 


Fig. ^our parameterizations of the NBB neutrino spectra by thick solid curves for several 
peak energies. For comparison, the corresponding original NBB spectra are shown by 
histograms. The parameterization reproduces well the main part of the NBB’s. Because 
it does not account for the secondary high-energy peak from decays (see 

Fig. i, we check that our main conclusions are not affected by those details (especially 
the background from CC events). 

We have not made parameterizations for the secondary (h^, z/g, z/e) beams. In the fol¬ 
lowing analysis we use the MC generated secondary beams at discrete energies (i^peak) 
and make interpolation for the needed i^peak values. Fluxes at different distances and for 
different species are obtained easily by multiplying the (2100 km/L)^ flux factor and the 
ratio of the cross section at a given E^^. 

4 Results 

In this section we present results of our numerical studies on physics potential of VLBL 
experiments by using the NBB’s from HIPA. 

First, we present our basic strategy of the analysis and explain our simplihed treatments 
of signals and backgrounds, and those of statistical and systematic errors. In the next 
subsection, we give our reference predictions for the results that may be obtained from the 
LBL experiment with HIPA and SK (L = 295km). In the latter two subsections, we give 
results of L = 2,100 km and L = 1, 200 km, respectively. 
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4.1 Signals, backgrounds and systematic errors 

In order to explore the physics potential of a VLBL experiment with HIPA at several base¬ 
line lengths, we make the following simple treatments in estimating the signals and the 
backgrounds of a future experiment. For a detector we envisage : 

• A 100 kton-level water-Cerenkov detector which has a capability of distinguish 
CC events from /i=*= CC events, but does not distinguish their charges. 


• We do not require capability of the detector to reconstruct the neutrino energy. 


Although water-Cerenkov detectors have the capability of measuring the three-momentum 
of the produced /i^ and as well as a part of hadronic activities, we do not make use 
of those information in our simplihed analysis. Instead we use only the total numbers of 
the produced /i^ and events from a NBB with a given peak energy. For each base-line 
length L, we study the impacts of splitting the assumed total experiment exposure of 1000 
kton-year (with 10^^ POT/year) into equi-partitioned runs of NBB’s at several peak en¬ 
ergies. We hnd that the use of two different-energy NBB’s improves the physics resolving 
power of the experiment signihcantly, but we have not found further improvements by 
splitting the experiment into more than two NBB’s. We therefore choose two appropriate 
NBB’s at each L, whose peak energies are chosen to make physics outputs (such as sensi¬ 
tivity to the neutrino mass hierarchy, sin^ and angles) signihcant. Optimum 

choice of NBB’s should depend on the model parameters 


5m 


2 

ATM 


5rrr 


sin 29 , sin 29^^-^ , sin 2,9^^^^^^ , 


). 


(4.1) 


especially on and sin^ which will be measured more accurately by K2K P], 

MINOS [0 and by HIPA-to-SK in the future. All our major hndings will not be 
affected by such details as long as appropriate NBB’s are chosen according to the data 
available at the time of the VLBL experiment. 

The signals of our analysis are the numbers of CC events and those of Ve CC events 
from the beam. They are calculated as 


/*-^peak 

A(/i, Epeak, L) = MNa / 

J 0 

/*£^peak 

N{e,Epeak,L) = MNa dEi,^{Epea.k) 

JO 






(4.2) 


(4.3) 


where the flux at a distance L is calculated from the parameterization at L = 2,100 km 


eq.(^) by multiplying the scale factor (2,100 km/L)^. The cross sections are obtained 
by assuming a pure water target. At low E^, the ratio of z/g and CC cross sections is 
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significant by different from unity, see eq.( p.4|) . Because of the vanishing of the NBB flux 
at low energies, our results are insensitive to the lower edge of the integration region. 
The probabilities and Py are calculated for the following model parameters ; 


20^^^ = 1.0 , 


,2 

sin^ 29= 0. 


'SOL 


5m^ = 3.5 X 10 eV , 

ATM ’ 

5m? = 10 X 10“^ eV^ , 


sin^ 29 ^ 


= 0.00, 0.02, 0.04, 0.06, 0.08, 0.10, 


^mns = 0°,90°,180°,270°, 
for the neutrino mass hierarchy I 

5ml^ = > 0, 5ml^ = > 0, (hierarchy I ) 

and for a constant matter density 

p = 3 g/cm^ 


(4.4a) 

(4.4b) 

(4.4c) 

(4.4d) 


(4.5) 


(4.6) 


throughout the base-line. We show in Fig. ^ the oscillation probabilities calculated for 
the above parameters, (at sir? 29^^^^^ = 0.1 and = 270° ) at three base-line lengths, 
L = 295 km (SK), L = 1, 200 km and L = 2 ,100 km. The NBB flux (xEy) chosen for our 
analysis are overlayed in each hgure. We use the result at SK (L = 295 km) expected for 
the low-energy NBB with (p^) = 2 GeV as a reference. 

The following background contributions to the ‘p’ and the ‘e’ events are accounted for : 

77(p, Lpeak, L)bg = ^(h; F^CC) -F iV(p; z/^CC : r ^ p), (4.7) 

7V(e, Lpeak, T)bg = N{e; u.CC) + iV(e; z/,CC : r ^ e) + iV(e; NC) 

-1- N{e; u^CC : t hadrons). (4.8) 

Here CC and Ug CC contributions are calculated by interpolating the numerical inte¬ 
grations 


iV(/; V CC) = MNa f dL,<F(_)(Lpeak)crf_‘(P(-) (_ 

J Vl VI Vl ^ Vl 


(4.9) 


for the discrete set of the MC simulations ET . Here ‘h(-)(Lpeak) and ^vXE^eak) stand for, 
respectively, the MC generated secondary and flux of the primarily beam. The 
survival probabilities are calculated for the same set of the model parameters, eg.(4.4). 
We hnd that contributions from oscillations from the background beams, such as W —> F^, 
are negligibly small and hence they are not counted. The contributions from r-lepton 
pure-leptonic decays are estimated as 

N{1- UrGC ■.T^l) = MNa J dEvHEpe.^)Pv^^vN?!'Br{T ^ , (4.10) 
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Figure 9: The neutrino oscillation probabilities at (a) L = 295 km, {h)L = 1, 200 km 
and (c) L = 2,100 km calculated for the parameter values ; sin^26*^^^ = 1.0, = 

3.5 X 10-W, sin2 203 oj^ = 0.8, 5ml^^ = 10 x lO-^eV^ = 0.1, 5^^^ = 270° 

and p = 3 g/cm^ with the hierarchy I. Overlayed are the NBB’s (flux times energy) used 


in our analysis. 
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where we adopt Br{T —> = 0.1737 and Br{T —>• ez/gZ/,-) = 0.1783 [^. The 10 


% errors in these branching fractions are accounted for as systematic errors. Because of 
the r-lepton threshold, r-backgrounds are signihcant only at high energies, NBB’s with 
h^peak ^ 4 GeV. Because they receive contribution from the small high-energy secondary 
peak due to Kaon decays, we use the interpolation of the results obtained for discrete set 
of MC generated fluxes. 

The ‘e’ events receive contributions from the NC events where produced vr'^’s mimic 
electron shower in the water-Cerenkov detector. By using the estimations from the SK 
experiments, we adopt 


iV(e;NC) = MNa J rf£„<*(£p..i.)<'=P./NC 


(4,11) 


with 


= 0.0055 ± 0.00055 


(4.12) 


The error in the above e/NC mis-identihcation probability is accounted for as a systematic 
error. The last term in eq. (|4.8|) accounts for the probability that the Uj. CC events with 
hadronic r-decays are counted as e-like events. In the absence of detailed study of such 
backgrounds, we use the same misidentideation probability for the NC events eq.( |4.12D and 
obtain 


N{e; z^T-CC : r —>• hadrons) = 


MNa / dE^^{Epe,.k)Pu, 




Brij —> z/T-hadrons)Pe/NC • (4.13) 


If the small misidentiheation probability of eq. (b:.12|) holds even for r —> hadron events, 
their background is only at the 2 % level of the r —>■ e background, iV(e; z^t-CC : r —*• e), 
and hence can safely be neglected. 

In addition, we account for the following two effects as the major part of the systematic 
uncertainty in the VLBL experiments. One is the uncertainty in the total flux of the 
neutrino beam, for which we adopt the estimate. 


/,„, = ! ±0.03, (4.14) 

common for all the high-energy NBB’s. We allocate an independent flux uncertainty of 3% 
for the low-energy NBB used for the SK experiment {L = 295 km) since it uses different 
optics. Finally, we allocate 3.3 % uncertainty in the matter density along the base-line. In 
our simplihed analysis, we use 

p = 3.0±0.1 g/cm^ (4.15) 

as a representative density and the uncertainty. 
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(a) L = 2 ,100 km 



^signal 

Abg 

iVtot 

-^^signal 

Mot 

beams 

n, GG 

NG 

NBB (Epeak = 6 GeV) 


202.5 

2.2 

15.2 

— 

219.9 

.92 

500kton-year 

N(e) 

126.5 

7.3 

15.9 

3.3 

153.0 

.83 

NBB (Epeak = 4 GeV) 

JV(^) 

612.5 

2.2 

3.5 

— 

618.2 

.99 

500kton-year 

N{e) 

66.4 

8.5 

3.7 

2.6 

81.2 

.82 


(b) L = 1, 200 km 



^signal 

Abg 

iVtot 

-^^signal 

Mot 

beams 

Ur GG 

NG 

NBB (Epeak = 5 GeV) 

N(p.) 

490.1 

2.2 

8.7 

— 

501.0 

.98 

500kton-year 

N{e) 

239.1 

10.5 

9.0 

9.4 

268.0 

.89 

NBB (Epeak = 3 GeV) 

NM 

413.8 

2.3 

0.0 

— 

416.1 

.99 

500kton-year 

N{e) 

186.1 

4.8 

0.0 

5.6 

196.5 

.95 


(c) L = 295 km 



^signal 

Nbg 

iVtot 

-^^signal 

Mot 

beams 

M GG 

NG 

NBB {{pr) = 2 GeV) 
lOOkton-year 


464.9 

7.8 

0.0 

— 

472.7 

.98 

N{e) 

161.3 

22.0 

0.0 

7.4 

190.7 

.85 


Table 6: Expected signals and backgronnds for /i-like and e-like events. The resnlts are 
shown for the parameters of eq.( |4.16| ) at = 270°, (a) 500kton-year at L = 2,100 km, 
(b) 500kton-year at L = 1, 200 km, (c) lOOkton-year at L = 295 km. 

Typical nnmbers of expected signals and backgronnds are tabnlated in Table 6 for the 
parameter set of eq. (|4.4|) at sin^ 26 '(.jjqq 2 = 0-1 (5^j^g = 270°. The nnmerical valnes are 

given for the following sets of experimental conditions ; 

(a) L = 2, 100 km 500 kton-year NBB(i7peak = 4 GeV), NBB(i7peak = 6 GeV) 

(b) L = 1,200 km 500 kton-year NBB(i7peak = 3 GeV), NBB(i7peak = 5 GeV) 

(c) L = 295 km 100 kton-year NBB((p 7 r) = 2 GeV) 

We note here that 100 kton-year at L = 295 km is what SK can gather in approximately 
5 years with 10^^ POT par year. 500 kton-year at longer distances can be accnmulated in 
5 years for a 100 kton detector with the same intensity beam. 

A few remarks are in order. A 100 kton-level detector at L = 2,100 km or 1,200 km 
can detect comparable nnmbers of /i-like and e-like events as SK (22.5 kton) at L = 295 
km. The backgronnds dne to secondary beams (mostly from the z/g beam ) are significant 
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in all cases for the e-like events. The NC background for the e-like events remain small at 


high energies if the estimate eq. (|4.12|) obtained from the K2K experiment at SK remain 
valid. We may expect gradual increase in the misidentification probability Pe/NC at high 
energies as the mean multiplicity of 7r° and charged particles grow. Finally the r-decay 
background can be significant only for NBB’s with Ppeak ^ 5 GeV. If a detector is capable 
of distinguishing r-events from the v^. and CC events, the overall fit quality improves 
slightly in the three-neutrino model, because the constraint -1- Pu^^ur — 1 

implies that the information obtained from the measurement always diminishes by 

contaminations from the r-events. As remarked in section 3.4, the CC event rates have 
been corrected for the high-energy secondary peak contributions by using the original MC 
generated NBB spectrum . 


4.2 Results for L = 295 km 


NBB <p„>=2Qey at L=295 km 



N{e)l^ OOkton-year 


Figure 10: CP phase dependence of N{e) and at SK(L = 295 km) for lOOkton-year 
with the NBB((p.n-) = 2 GeV). = 0° (solid-circle), 90° (solid-square), 180° (open- 
circle), and 270° (open-square). The results are shown for sin^ = 1.0, = 

3.5 X 10“^eV^, sin^ 26*gQL = 0.8, = 10 x 10“®eV^, sin^ 26 *(.jjqq 2 = 0.1, and p = 3 

g/cm^. The predictions for the four types of the neutrino mass hierarchies (Fig. |l|) are 
depicted as I, II, III and IV. 


In Fig. 1^ we show the number of CC events, N{p), and that of Ue CC events, 
N{e), expected at SK (L = 295 km) for 100 kton-year. The NBB with {p^) = 2 GeV, 
NBB((pjr) = 2 GeV), with lO^^POT per year has been assumed for simplicity. The three- 


neutrino model parameters of eq. (|4.4|) at sin^ 29^ 


= 0.1 are assumed, and we take 
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p = 3g/cm^, eq. (|4.6|) . The predictions then depend only on the CP phase parameter 
and when we vary from 0° to 360°, we have a circle on the plane of N{p) vs N[e). 
The four representative cases, = 0°, 90°, 180° and 270°, are marked by solid-circle, 
solid-square, open-circle, and open-square, respectively. The four possible mass hierarchy 
cases of Fig. | are depicted as I, II, III and IV. 

Only the expected event-numbers from the beam, eq. ([4.2|) and eq. (|4.3|) , are counted 
in Fig. 0 so that we can read off the ultimate sensitivity of the experiment from the 
hgure. Statistical errors for such an experiment are shown for the = 270° point in the 
mass hierarchy I circle. We can learn from the hgure that if we know all the parameters 
of the three-neutrino model except for the mass hierarchy and if we know the neutrino- 
beam hux exactly, then there is a possibility of distinguishing the mass hierarchy I from 
III. In practice, the LBL experiment between HIPA and SK can constrain mainly 
and sin^26'^pj^ from N{p), and from N{e) [^. Fig. shows that those 

measurements should suffer from uncertainties in the remaining parameters of the three- 
neutrino model, the neutrino mass hierarchy cases and explicitly shown, as 

well as on the solar-neutrino oscillation parameters, 5rn? and sin^ 26'gQp. Next generation 


of the solar-neutrino observation experiments and KamLAND experiment [Q may 
further constrain the latter two parameters, but the mass-hierarchy (between I and III) and 
should be determined by the next generation of accelerator-based LBL experiments. 
It is hence necessary that all the results from the LBL experiment between HIPA and SK 
should be expressed as constraints on the three primary parameters, , sin^ ‘^^atm > 


and sin^ 26 '(^jjqq 2 5 which depend slightly on the three remaining parameters of the three- 
neutrino model, 5^gQL, sin^ 26^3^^ and mass-hierarchy cases. In the 

following subsections, we show that the data obtained from the LBL experiment between 
HIPA and SK {L = 295 km) are useful in determining the neutrino mass hierarchy, and in 
some cases even when they are combined with the data from a VLBL experiments 
(L = 2,100 km or 1,200 km) with higher-energy neutrino beams from HIPA. 

Before moving on to studying physics potential of VLBL experiments, it is worth noting 
that the predictions for the mass hierarchy IV in Fig. 0 represents the prediction for the 




oscillation probabilities in the hierarchy I, according to the theorem eq.( |2.57|) , 
once the scales are connected by the factor aiyiN)/aiviN). By comparing the depen¬ 
dences of the circle I and circle IV, we can clearly see that and interchange 

approximately by exchanging 5^j^g = 90° and 270°. The comparison of and oscilla¬ 
tion experiments at around L = 295 km hence has a potential of discovering CP-violation 
in the lepton sector. However, determination of the angle by using the beam from 


HIPA needs a much bigger detector than SK |19, 46 
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4.3 Results for L = 2,100 km 


Fig. 0 shows the numbers of CC events and that of z/g CC events expected at the 
base-line length of L = 2,100 km from HIPA with 500 kton-year. The expected signal 
event numbers are shown for (a) the NBB with F^peak = 4 GeV and for (b) the NBB with 
F'peak = 6 GeV. The parameters of the three neutrino model and the matter density are 
taken exactly the same as in Fig. 0: 


sin^ 29^,^ = 1.0 , = 3.5 x lO'^ eV^ , 

sin^ = 0.8 , = 10 x 10“^ eV^ , 

SOL ’ SOL ’ 

sin^20cHooz = 0-G 


(4.16a) 

(4.16b) 

(4.16c) 

(4.16d) 

(4.16e) 


ATM 


ATM 


CHOOZ 



p = 3 g/cw? . 


The predictions for the four neutrino mass hierarchy cases (Fig. |l|) are shown by separate 


circles when the GP-phase angle is allowed to vary freely. The predictions for the 
four representative phase values are shown by solid-circle (^^^g = 0°), solid-square (90°), 
open-circle (180°), and open-square (270°). 


(a) NBB (Epeak=4GeV) at L=2100 km (b) NBB (Epeak=6GeV) at L=2100 km 
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Figure 11: GP phase dependence of N{e) and N{p) at L = 2100 km for 500 kton-year with 
(a) = 4 GeV) and (b) NBB(Bp„^ = 6 GeV), = 0» (solid-circle), 90» (solid- 

square), 180° (open-circle), and 270° (open-square). The results are shown for 5m? = 


3.5 X 10“^eV^, sin^ = 10 ^ 10“^eV^, sin^ = 0.8, sin^ ^ 0.1, 

and p = 3 g/crv?. The predictions for the four types of the neutrino mass hierarchies (Fig. 
are depicted as I, II, III and IV. 
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When comparing with the L = 295 km case (Fig. 0 ). it is most striking to find that the 
predictions for the CC events, N (e), differs by a factor of 5 or even larger in magnitude 
between the neutrino mass hierarchy I and III. This is because of the enhancement of the 
matter effect at high energies as depicted in Fig. This striking sensitivity of the prob¬ 
ability on the mass hierarchy cases is the bases of the capability of distinguishing 

the cases in VLBL experiments by using the HIPA beam. On the other hand, we will find 
that the 5 % level differences in N{^) between the mass hierarchy cases are not useful for 
this purpose because depends strongly on the parameters Sml and sin^ 


and also because A^(/i) suffers from the neutrino-beam flux uncertainty eq. ([4.14|) . 

We have examined NBB’s with various peak energies, and find that the NBB with 
-Fpeak = 6 GeV (see Fig. |^) makes N{e) largest while keeping N{fi) small for the atmospheric- 
neutrino oscillation parameter of eg. ( [4. 16^) . The ration N{e)/N{fi) can be as large as 1/2 
if sin^ = 0.1 (eq. (|4.16d) ) for the hierarchy case I. The NBB with F^peak = 4 GeV is 

then chosen because it has the dependence of N{e) which is ‘orthogonal’ to that of 
the i^peak = 6 GeV case. By comparing the hierarchy I circles of Fig. |TT](a) and (b), we find 
that for the NBB with Epe^k = 4 GeV, N{e) is largest at around = 180° (open-circle) 
and smallest at around = 0° (solid-circle), whereas for the NBB with i^peak = 6 GeV, 
N{e) is largest at around 5^j^g = 270° (open-square) and smallest at around hj^j^g = 90° 
(solid-square). 

In order to understand the dependences of V(p), N{e) on i^peak (the peak energy of the 
NBB) and on and in particular to probe the reason for this peculiar dependence 
of the expected number of GG events on F^peak, we show in Fig. the oscillation 
probabilities and Pu^^ue plotted against the neutrino energy, E,^. The three plots 

in the left-hand-side are for and the three in the right-hand-side are for Py^^y^. In 

all the plots, the predictions of the three-neutrino model with the parameters of eq. (|4.16|) 
are shown by thick lines for the neutrino mass hierarchy I. In the upper, middle and bottom 
figures, the predictions for the same parameters are shown by thin lines for the hierarchy 
II, III, and IV, respectively. In all the plots the predictions for h^^^g = 0°, 90°, 180°, 270° 
are shown by solid, long-dashed, short-dashed, and dot-dashed lines, respectively. It is 
clear from these figures that the transition probability Py^^y^ is more sensitive than the 
reduction probability Py^^y^ on the neutrino mass hierarchy and ^^ns- chosen 

parameters of eq.( ^4.16| ), Py^^y^ hits zero at around Ey = 6 GeV at F = 2,100 km, where 
the transition probability Py^^y^ becomes largest in the interval 4 GeV < Ey < Q GeV. 
Because the NBB with i^peak = 6 GeV is largest at the peak and has significant tail down 
to around F^i. = 4 GeV, see Fig. ||(c) for the NBB shape, we find largest N{e) and smallest 
N{fi) for this NBB. At the same time we clearly see that in the interval 4 GeV < Ey < 6 
GeV, Py^^y^ is largest at = 270° (dot-dashed line), and smallest at 5^^g = 90° (long- 
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dashed line). Just below Ey = A GeV the ordering between the = 270° (dot-dashed) 
curve and the = 180° (short-dashed) curves and that between the J^j^g = 0° (solid) 

curve and the = 90° (long-dashed) curve is reversed. Consequently the NBB with 

-E'peak = 4 GeV gives largest iV(e) at around J^^^g = 180° and smallest iV(e) at around 
Jmns = 0°) among the four representative angles. The overall normalization of N{e) is 
smaller, and that of N (/i) is larger for NBB (T^peak = 4 GeV) than the predictions of NBB 
(Epeak = 6 GeV). 

In Fig. 0, we show the statistical errors of the iV(/r) and N{e) measurements at 500 
kton-year on the J^^^g = 270° point for the hierarchy case I. The size of the error bars 
suggest that a 100 kton-level detector is needed to explore the model parameters in VLBL 
experiments at L ~ 2,100 km. It also tells us that such detector has the potential of 
discriminating the neutrino mass hierarchies and constraining the angle in a certain 
region of the three-neutrino model parameter space. A more careful error analysis that 
accounts for backgrounds and systematic errors is given in the next subsection. 

We are now ready to study the physics capability of such VLBL experiments in some 
detail. In Fig. |T^, we show the expected numbers of signal events, V(/i) and V(e), for the 
same NBB’s and for the same volume of 500 kton-year as in Fig. 0 The circles of Fig. 0 
are now made of 36 J^^^g points, J^^^g = n x 10° for n = 1 to 36, for each set of the model 
parameters. The common parameters for all the points are 


‘^^ATM = 1-0 ) = 3-5 X 10 ^ eV^ , (4.17a) 

sin^ = 0.00, 0.02, 0.04, 0.06, 0.08, 0.10, (4.17b) 

Jmns = n X 10° (n = 1 to 36), (4.17c) 

p = 3 g/cm^ . (4.17d) 


The predictions for the six sin^ 29^ 


cases eq.(|4.17bD can be recognized as hve circles 


and one point with decreasing N{e) values, as 8111^29^^^^^^ decreases from 0.1 to 0.0. The 
remaining two parameters are constrained by the solar-neutrino oscillation experiments, 
and we chose the following representative parameter sets for the three possible solutions 
to the solar-neutrino dehcit anomaly : 


LMA : 

sin^ 20go^ 

o 

bo 


f 15 X 10-^ eV^ , 

“ j 5 X 10-5 eV^ , 

(4.18a) 

SMA : 

sin^ 20go^ 

= 7 X 10-^ 


= 5 X 10-5 eV^ , 

(4.18b) 

VO : 

sin^ 20go^ 

= 0.7, 


= 7 X 10-1^ eV^ . 

(4.18c) 


The predictions of the four types of the neutrino mass hierarchy, I, II, III, IV are indicated 
explicitly. 
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(a) NBB( Epg3^=4GeV ) at L=2100 km 


(b) NBB( 6GeV ) at L=2100 km 
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Figure 13: The neutrino parameter dependences of the expected numbers of z/g CC and 
CC events, N{e) and iV(/i), respectively, for the NBB with -Epgak = 4 GeV (a) and 6 GeV (b) 
with 500 kton-year at L = 2,100 km. The predictions are shown for the LMA, SMA, and 
VO scenarios of the solar neutrino oscillations and for the four neutrino mass hierarchies I 
to IV. All the numbers are calculated for = 3.5 x 10“^eV^ and sin^ 29 = 1.0 at 

ATM ATM 

= 0.02 X m (m = 0 to 5) and = 10° x n (n = 1 to 36). The five larger 
circles in each hierarchy show the predictions of the LMA with = 15 x 10“®eV^ 

and the hve smaller circles are for the LMA with = 5 x 10“®eV^, both for m = 1 

to 5. dependence is not recognized for the predictions of SMA and VO, and also at 
sin^ = 0 (m = 0). The VO predictions are the same for the hierarchy I and II, and 

III and IV. 


The solid-circle points in Fig. show the predictions of all the models with neutrino 
mass hierarchy 1. They reside in the corner at large N{e) and small The five 

grand circles with smallest A(/i) give the predictions of the LMA solution for = 

15 X 10“® eV^, and those with larger V(p) are for the LMA with = 5 x 10“® 

eV^. It is clearly seen that the dependence (the size of the grand circles) is larger 
for larger and for larger sin^ ‘^^chooz expected. It is worth pointing out that 

the LMA scenario with = 15 x 10“® eV^ predicts non-zero N{e),N{e) ~ 5/500 

kton-year even when sin^ 26 *g,jjQQ 2 = 0 - This is because the ‘higher’ oscillation modes of the 
three-neutrino model grow as ^is^s. The predictions of the SMA parameters 

appear just above the upper LMA grand circles, where the h^^^g dependence (the size of 
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the grand circle) diminishes to zero for each siia^ value of eq. (^4.17b|) . Just hve 

solid-circle points appear in Fig. (a) and (b) for the SMA parameters because the 
Jmns” dependence diminishes as sin^ 203 ^ 1 ^ are both signihcantly smaller 

than unity. This is expected from the Hamiltonian eq.( p.55a| ) that governs the neutrino 
oscillation in matter, because the -dependence of observables diminished whenever 
the second terms proportional to 5m\2 /much smaller than the hrst term. The 
predictions of the VO parameters, which are given by the solid-diamonds, are shown just 
above those of the SMA parameters. Here the J^^^g dependence vanishes because of the 
extreme smallness of 5m? /dm? . With the same token, we cannot distinguish the VO 
predictions of the neutrino mass hierarchy I and H. Because the magnitude of is 

so small, the sign of 5m\2 = does not have observable consequences in terrestrial 

LBL experiments. 

The predictions of the four scenarios of the solar neutrino oscillations (VO, SMA, 
and two cases of LMA in eq. (|4.18|) ) with the neutrino mass hierarchy H (Jm ^2 = 

—5m‘^^^,5mi^ = are shown by open-circle points, which are located in the corner 

of large N{e) large As explained above, the predictions of the VO scenario do not 

differ from those for the hierarchy I. The predictions of the SMA scenario appear slightly 
above those of VO. The two LMA scenarios predict larger V(/i) and give visible J^j^g de¬ 
pendences, which lead to hve grand circles each for the hve discrete values of assumed 
sin^ Again the -dependence (the size of the grand circles) is larger for the 

large sin^ 26 *^,jjqq 2 large 5m'^^^. Summarizing the predictions of the four scenarios in 


the neutrino mass hierarchies I and H, both of which have Jm?, = 6m‘^ > 0, N(e) 

grows linearly with increasing sin^ ‘^^caooz scenarios and for both hierarchies. 

The predictions of the hierarchy H diher from those of the hierarchy I only by slightly 
larger N (/i) in each scenarios. The difference in N[fi) is largest in the LMA scenarios with 
5m‘^o^ = 15 X 10“® eV^, for which the hierarchy H predicts about 30 % (10 %) larger V(p) 
than the predictions of the hierarchy I for the NBB with i^peak = 6 GeV (4 GeV). 

All the predictions of the hierarchy cases HI and IV have very small N{e), as we ex¬ 
plained in section 2.4 for the hierarchy HI. This is because the two mass hierarchies have 
common. The predictions of the hierarchy HI are shown by solid-squares 
and those of the hierarchy IV are shown by open-squares. As in the case of hierarchy I 
vs H, the predictions of the VO scenario do not show visible dependence on the sign of 
5ml2 = common predictions for the hierarchy HI and IV are shown by 

open-diamonds. Even though compressed to the very small N{e) region, the largest N{e) 


point is for sin^ 26', 


= 0.1. The predictions of the SMA scenario appear just above 


(below) the VO points for the hierarchy HI (IV). In the LMA scenario, the predictions for 
V(/i) increases as grows for the hierarchy HI, while N{fi) decreases with for 
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the hierarchy IV. By noting that the > z/g oscillation probability of the hierarchy IV is 
approximately that of the > Pg oscillation probability of the hierarchy I, see eq. (|2.57|) , 
and by noting that the Pg CC cross section is a factor of three smaller than the CC 
cross section, we can conclnde from the fignres that the P^ —>■ Pg oscillation experiments 
by using the P^ beam from HIPA is not effective if the neutrino mass hierarchy is indeed 


type I [33 


4.4 A case study of semi-quantitative analysis 


Now that the intrinsic sensitivity of the very simple observables iV(/i, ii^peak) and N{e, i^peak) 
on the three-neutrino model parameters, sin^ and the neutrino mass 

hierarchy cases are shown clearly in Fig. 0 . we would like to examine the capability of 
such VLBL experiments in determining the model parameters. The following 4 questions 
are of our concern: 


1 . 

2 . 

3. 

4. 


Can we distinguish the neutrino mass hierarchy cases ? 

Can we distinguish the solar-neutrino oscillation scenarios 

Can we measure the two unknown parameters of the model, sin^ ‘^(^caooz 5^^g ? 
How much can we improve the measurements of 6ml and sin^ ‘^^atm 


In this subsection we examine these questions for a VLBL experiment at L = 2,100 km in 
combination with the LBL experiment of L = 295 km between HIPA and SK. 


For definiteness, we assume that the VLBL experiment at L = 2,100 km accumulates 
500 kton-year each with the NBB (i^peak = 6 GeV) and NBB (F^peak = 4 GeV). As for the 
LBL experiment at L = 295 km, we assume that 100 kton-year data is obtained for the 
NBB ((Ptt) = 2 GeV). Although this latter assumption does not agree with the present 
plan of the HIPA-to-SK project where the GAB may be used in the first stage, we find 
that essentially the same final results follow as long as the total amount of data at L = 
295 km is of the order of 500 kton-year. The expected numbers of p-like and e-like events 


V(/i) — V(p, F'peak, L) + N{fi, i^peak, F)bG , (4.19) 

N{e) = N (e, F'peak, L) + N{e, F^peak, F)bg , (4.20) 


are tabulated in Table ^ for the parameters of eq. (j4.16| ) at = 0°. We see from 

the Table that roughly the same number of > z/g signal events is expected for an 
experiment at F = 2,100 km with 500 kton-year and for an experiment at F = 295 km 
with 100 kton-year. This agrees with the naive scaling low of V ~ 1/F at same LjEy. 

Our program proceeds as follows. For a given set of the three-neutrino model pa¬ 
rameters, we calculate predictions of V(/i) and N{e), including both the signal and the 
background, with a parameterized neutrino-beam flux and a constant matter density of 
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the earth (p = 3g/cm^). We assume that the detection efficiencies of the p-like and e-like 
events are 100% for simplicity. The statistical errors of each predictions are then simply 
the square roots of the observed numbers of events, N{fi) and N{e). In our analysis, we 
account for the following systematic errors: 


• Overall flux normalization error of 3% 

• Uncertainty in the matter density along the baseline of 3.3% 

• Relative uncertainty in the misidentihcation probability Pe/NC of 10% 

• Relative uncertainty in the r-backgrounds of 10% 


The function of the fit to the LBL experiments may then be expressed as the sum 


= 2,100 km) + x^{L = 295 km), 


(4.21) 


where the hrst term is 


X^{L = 2,100 km) 


) I J FLUX 

£^peak, kton-year 


/flux-^^**(/^) - V , - iV*’--(e) V' 


a(p) 


a{e) 


7flux-1^ ^fP-3 


0.03 


0.1 


(4.22) 


o-(/i) = ^7Vi™®(/i) -h (0.1iV*™'=(/i; UrCC)f , (4.23) 

a(e) = i/w™®(e) (0.1iV*™'=(e; NC))^ -h (0.1iV*™®(e; z/^CC))7 (4.24) 

Here the is a function of the parameters of the three-neutrino model (the two mass- 
squared differences, three angles and one phase), the flux normalization factor /p^ux 
matter density p. We assign the overall flux normalization error of 3% which is common 
for all high-energy NBB(Ppeak)’s, including the secondary beams. Besides this common 
flux normalization error, the individual error of N{p) is the sum of the statistical error 

and the systematic error coming from the uncertainty in the r-background. The error of 

7V(e) is a sum of the statistical error and the systematic errors coming from the uncertainty 
in the e/NC misidentihcation probability and the r-background. The function for the 
HIPA-to-SK experiment is obtained similarly where the 3% hux normalization error and 
the 3.3% matter density error are assumed to be independent from those for the L = 2, 100 
km experiment. 
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(a) Hierarchy I (true) vs III (fit) (b) Hierarchy I (true) vs III (fit) 



Figure 14: Minimum as functions of the fitting parameter sin^ 26*(3 .jjqq 2 by assuming the 
hierarchy III, when the mean values of iV(/i, -Speak) and ./V(e, -Speak) are calculated for the 
LMA points with hierarchy I (eg. ( [4.2(j| )). In total 12 cases of the input data are labeled 
by the input (‘true’) values of sin^ = 0.02 (thin lines), 0.04 (medium-thick lines), 

0.1 (thick lines) and = 0° (solid lines), 90° (long-dashed lines), 180° (short-dashed 
lines), 270° (dot-dashed lines). The minimum of the function is found by assuming the 
hierarchy III for a since sin^ value, by varying the parameters and sin^ 

within the LMA allowed region, and the remaining three parameters, sin^ 26^^^ 

and freely (eq. (|4.27|) ). (a) Result with 500 kton-year each for NBB(.Speak = 6 GeV) 
and NBB(.Speak = 4 GeV) at S = 2,100 km. (b) In addition, 100 kton-year data from 
NBB((p^) = 2 GeV) at S = 295 km are included in the £t. 

We hrst examine the capability of the VLBL experiment in distinguishing the neutrino 
mass hierarchy cases I and III. For this purpose, we calculate the expected numbers of 
signals and backgrounds for a certain set of the three-neutrino model parameters with the 
hierarchy I, and the examine if the generated data can be interpreted under the assumption 
of the hierarchy III. 

4.4.1 Neutrino mass hierarchy 

We show in Fig. |T^ the minimum as functions of the parameter sin^ ‘^^caooz by assuming 
the hierarchy III, when the mean values of N{fi, Epeak) and N{e, -Speak) are calculated for 
the LMA points with the hierarchy 1. The results in Fig. 0 (a) and (b) are given for the 
following sets of experimental conditions (A) and (B), respectively: 

(A) 500 kton-year each for NBB (.Speak = 6 GeV) and NBB (.Speak = 4 GeV) [a ok \ 

ofr_oinnUw. (4.25aj 
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(4.25b) 


(B) In addition to (A), 100 kton-year data from NBB((p 7 r) = 2 GeV) 
at L = 295 km are included in the fit. 


We take the 12 sets of A^(/r, ii^peak) and A^(e, ii^peak), which are generated for the LMA 
parameters of, eq.(^^) at = 0.02, 0.04, 0.1, with the matter density p = 3 

g/cm^ by assuming the hierarchy I. Summing up, the expected numbers of signals and 
background events are calculated for the following input (‘true’) parameters: 


in2 20true 
ATM 

= 1 . 0 , 

trne ^ 3 5 ^ ^q-3 y 2 

ATM ’ 

(4.26a) 

sin 2 2 C/ 

o 

bo 

^ ^q-S gy 2 ^ 

(4.26b) 


= 0.02, 0.04, 0.10, 

Cns = 0°’ 90°, 180°, 270° , 

(4.26c) 

5ml^ 

= > 

, (hierarchy I) 

(4.26d) 

£true 

''FLUX 

= 1 . 0 , 

^true ^ 3g/cjn3 _ 

(4.26e) 


The 12 cases of the input data sets are labeled by the input (‘true’) values of sin^26*^™®Q^ 
= 0.02 (thin lines), 0.04 (medium-thick lines), 0.1 (thick lines) and 5^“® = 0° (solid lines), 
90'’ (long-dashed lines), 180'’ (short-dashed lines), 270'’ (dot-dashed lines). 

The fit hs-s been performed by assuming the hierarchy III. The minimum of the 
function is found for a sin^ ‘^^chooz value in the range below 0 . 12 , by varying the parameters 
hrUgoL cind sin^26*gQL within the LMA allowed region, eq. (|2.32|) , and the remaining three 
parameters, sin^ 20^.^,^ and freely. The uncertainties in the total fluxes of 

the neutrino beams, eq. (|4.14|) , and the matter density, eq. (|4.15|) , are taken into account. 
Summing up, the fitting parameters used to obtain the functions of Fig. IT3 are : 


: free, ’ (4.27a) 

sin^ 20go^ = 0.7 - 0.9 , 6ml^^ = (3 - 15) x 10“^ eV^ , (LMA) (4.27b) 

sin 20 (^jjqq 2 • free, ^mns • fr^e, (4.27c) 

5ml2 = Sml^^ , dmls = , (hierarchy III) (4.27d) 

/flux : free, p : free. (4.27e) 


In Fig. 1^ (a), we show the resulting Xmin froni the VLBL experiments at L = 2,100 km. 
The minimum x^ for = 0.1 is always larger than 30 for sin^ 26 *(.jjqq2 < 0.12. 

We can hence distinguish the neutrino mass hierarchy I from III. It is also possible to 


make the distinction at more than 3a level for the cases with = 180'’ and 270° at 

MNS 

sin^26**™® , = 0.04 whereas the distinction disappears for = 0° and 90°. This is 
because, the hierarchy I predictions of A^(e)’s for = 0° and 90° at sin^ ‘^^chooz ~ 
can be reproduced by the hierarchy-III model if we choose larger sin^ 26 *(.jjqq 2 ( ^ 0.08); see 


Fig. 11 and Fig. 13 


42 












In Fig. |TJ(b), the minimum values are shown when data from HIPA-to-SK experiment 
are added. The remarkable difference between Fig. 0 ( a) and (b) is found for the small 
sin^ cases when the htting parameter sin^ ‘^(^caooz large. This can be explained 

as follows. In the HIPA-to-SK LBL experiment, the predicted A^(e, i^peak) of the hierarchy 
III is not much smaller than that of the hierarchy I (see Fig. |^. In particular, N{e, T^peak) 
calculated for the hierarchy I at sin^ ‘^^cnooz ^ signihcantly smaller than that for 

the hierarchy III at sin^ ‘^^caooz >0.1. This leads to the enhancement of the minimum 
in Fig. |T^(b) at large sin^ We hnd that the data obtained from the HIPA-to-SK 

experiments are useful, which allow us to determine the neutrino mass hierarchy (between 
I and HI) for all four values of at 3a level when sin^ ‘^^ckooz ^0.04, and at one-sigma 
level when sin^ ^ 0.02. 

Let us now study the possibility of distinguishing the solar-neutrino oscillation scenarios 
by using the VLBL experiments with HIPA. Because the predictions of the SMA and VO 
scenarios differ very little in Fig. 1^ we compare the predictions of the SMA and LMA 


scenarios. 


(a) LMA (true) vs SMA (fit) 


(b) LMA (true) vs SMA (fit) 




Figure 15: Minimum as functions of the parameter sin^ ‘^^caooz- "^1^® mean values of the 
input data, V(/i, T^peak) and N{e, F^peak), are calculated for the 8 LMA points: sin^ ‘^^tnooz 
= 0.06 (thin lines), 0.1 (thick lines) and = 0° (solid lines), 90° (long-dashed lines), 
180° (short-dashed lines), 270° (dot-dashed lines). The £t has been performed by 
assuming the SMA scenario for and sin^ ‘^^sol varying the remaining four pa¬ 
rameters, sin^ 20^.^.,^, and sin^ 20^^^^^. (a) With 500 kton-year data each with 

NBB(i7peak = 6 GeV) and NBB(i7peak = 4 GeV) at L = 2,100 km. (b) In addition, 100 
kton-year data from NBB((p 7 r) = 2 GeV) at L = 295 km are included. 
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4.4.2 LMA v.s. SMA/LOW/VO 


We show in Fig. the minimum as functions of the parameter sin^ ‘^^caooz t’Y assuming 
the SMA scenario, when the mean values of A^(/x, F'peak) and A^(e, F^peak) are calculated for 
the LMA point with hierarchy I (eq. ([4.26|) ). The results of Fig. (a) are obtained from 
the L = 2 ,100 km VLBL experiments only, (b) are obtained after adding the HIPA-to-SK 
data; see (A) and (B) in eq. (|4.25|) . 

We use the same input parameters as in eq. (|4.26|) , but only two cases of sin^ ‘^^cnooz ’ 
0.06 and 0.1, are examined. The 8 cases of the input data are again labeled by the input 
(true) values of sin^26*^™®^^ = 0.06 (thin lines), 0.1 (thick lines) and = 0° (solid lines), 
90° (long-dashed lines), 180° (short-dashed lines), 270° (dot-dashed lines). 

The fif has been performed by assuming the SMA solution, eq. (|4.18b|) , with the 
hierarchy I. The other fitting parameters in the three-neutrino model are taken for the 
same as in Fig. B Summing up, the parameters used in the fit are 


sin" 20^^^ 

: free, 

ATM 

: free, 

(4.28a) 

sin" 

= 7 X 10■^ 


= 5 X 10"® eV^, (SMA) 

(4.28b) 

sin 20 ^^jjqp 2 

: free, 

A 

'^MNS 

: free, 

(4.28c) 


= ’ 

: 

= ’ (hierarchy I) 

(4.28d) 

/FLUX 

: free, 

p : free. 

(4.28e) 


In Fig. |^(a), the minimum becomes almost zero when = 90° or 270°, whereas 
the results with the 5*™® = 0° and 180° give non-zero values of the minimum v^. This is 
because the two beams, NBB (i7peak = 4 GeV) and NBB(£'peak = 6 GeV), give ‘orthogonal’ 
(5 mns dependences of N(e, i7peak) in Fig. |T^ For = 90° and 270°, the predictions of the 
LMA scenario for A(e,4 GeV) and N(e,6 GeV) can be reproduced in the SMA scenarios 
where the fitting parameter sin^ is approximately chosen, but such reproduction 

does not work perfectly for the = 0° and 180°. Although the minimum x^ value is 
raised in Fig. |T^(b) after combining the HIPA-to-SK data, it seems difficult to distinguish 
the LMA scenario from the SMA scenario. 

The x^ fit by assuming the LOW or VO scenarios give the same results as those obtained 


by assuming the SMA scenario. This is understood from the Hamiltonian eq. (|2.55a|) that 
governs the neutrino oscillation in the vacuum. The effects of solar-neutrino oscillation 
modes are proportional to the term with the LOW and VO scenarios it is 

the extreme smallness of /m^ that makes the effects unobservable. In case of the 

SMA scenario, the effects are also tiny because the mixing matrix elements multiplying 
the Im? factor are also very small. 

SOL/ ATM 
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Figure 16: Regions allowed by the VLBL experiment at L =2,100 km are shown in the 
plain of sin^ and The experimental conditions are given in eq. (|4.25aj) . The 

input data are calculated for the LMA parameters, eg. ( |4-26| ), at sin^ 20^™®^^ = 0.1 and 
for four values of 5^“®, 0°, 90°, 180°, and 270°. In each figure, the input parameter point 
(sin^ is shown by a solid-circle, and the regions where x^min <1)4, and 9 are 

depicted by solid, dashed, and dotted boundaries, respectively. 


4.4.3 


sim 26 and 5 

CHOOZ MNS 


Fig. 1^ shows the regions in the sin^ 26 *(.jjqq2 ^mns plane which are allowed by the 
VLBL experiment at L =2,100 km with 500 kton-year each for NBB(i7peak = 6GeV) and 
NBB(i?peak = 4GeV); eq.( [4.25a|) . The input data are calculated for the LMA parameters, 
eq. ([4.26| ), at sin^ 26*^™®^^ = 0.1 and for four values of , 0°, 90°, 180°, and 270°. The 
fit has been performed by assuming that and sin^ are in the LMA region 

eg.( 2.32 ) while the rest of the parameters are freely varied; 


: free, ’ (4.29a) 

sin^ 20goL = 0.7 - 0.9 , = (3 - 15) x lO"® eV^ , (LMA) (4.29b) 
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Figure 17: Regions allowed by the VLBL experiment at L =2,100 km and the LBL 
experiment at L =295 km, with the experimental condition of eq. (|4.25b|) . The symbols are 
the same as those in Fig. ITO. 


1 

CHOOZ 

: free, 

<5mns : free. 

(4.29c) 

6ml2 

= ^"'SOL > 

6ml^ = , (hierarchy I) 

(4.29d) 

f FLUX 

: free, 

p : free. 

(4.29e) 


In each hgure, the input parameter point (sin^ ‘^^ckooz ’ ^mns) shown by a solid-circle, and 
the regions where Xmin <1)4, and 9 are depicted by solid, dashed, and dotted boundaries, 
respectively. 

A few comments are in order. From the top-right hgure for = 90°, we learn that 
is not constrained by this experiment. The reason can be qualitatively understood by 
studying the LMA predictions shown in Fig. 0 and Fig. |T^. Note, however, that the LMA 
cases shown in Fig. |^are for = 5x 10“^eV^ and 15x 10“^eV^, while the input data in 

our example are obtained for = 10 x 10“®eV^. An appropriate interpolation between 

the two cases is hence needed. From Fig. o we hnd that = 90° (solid-square) points 
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Figure 18: The same as Fig. Q. but when the input data are calculated for sin^ 20 
0.06. 


true 

CHOOZ 


for the hierarchy I he in the lower N (e) corner of the grand circle made of the predictions of 
all (5 mns- On the other hand, we can tell from Fig. that the same values of low N{e) can 
be obtained for the other 6, 


MNS 


values typically by reducing the sin^ ‘^^caooz Parameter. In 
case of = 270°, both and sin^ are constrained in the bottom-right hgure, 

MWo iViiNO OFIOOZ 

because = 270° gives the largest N{e) along the relevant grand circle in Fig. pT|(b). 
When are chosen differently, we need to increase sin^ “^^cftonz niake N{e) large, but 


it then gets into conflict with the CHOOZ bound eq.( p.27D . 

The left hgures of Fig. p!8|show that can be constrained at la level when 5^“® = 0° 
(top-left), or at 2a level when 180° (bottom-left). The reason for this behavior is more 
subtle. We hnd that it is essentially the ratio of A^(e, 4GeV)/iV(e, 6GeV) which is difficult 
to reproduce with the other In Fig. 0, we can see among the four representative 

^MNs predictions that iV(e, 4GeV)/iV(e, 6GeV) is smallest at = 0° (solid-circle) and it 
is largest at = 180° (open-circle), while this ratio is almost the same for = 90° 
(solid-square) and 5^j^g = 270° (open-square). Because it is the energy-dependence of N{e) 
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^true _ no 

^MNS “ 

O 

o 

180° 

270° 

(xlO-3) (eV^) 

(A) 

'J-'-'U-O.lO 

3.5071):)^ 


3.507D:1^ 

(B) 

“o-EnTTrUT 

'J-"-'U-0.06 

3.5071)7? 

“o-EpTonir 

3.507D:1? 

h (sin220^^^) 

(A) 

> 0.976 

> 0.976 

> 0.975 

> 0.976 

(B) 

> 0.991 

> 0.991 

> 0.991 

> 0.991 


Table 7: The mean and the one-sigma errors of and the one-sigma lower bound 

of sin^ 20^*^ when the input data are calculated for the four LMA points (eq. (^4.26|) at 
sin^ = 0.06) and the fits are performed by assuming the LMA scenarios. (A) 

Results with 500 kton-year each for NBB(i7peak = 6 GeV) and NBB(i7peak = 4 GeV) at 
L = 2,100 km. (B) In addition, 100 kton-year data from NBB((p 7 r) = 2 GeV) at L = 295 
km are included in the £t. 


that has some discriminating power for detectors with the capability of measuring 


neutrino energy 1^ may have better sensitivity for the 6 angle. 


In Fig. 1^ we show the same plots as Fig. where in addition to the data from the 
VLBL experiment at L =2,100 km, we also include the 100 kton-year data from the L =295 
km experiment; see eq. ([4.25b|) . Although the area of Xmin 4, and 9 regions decrease 
signihcantly, the qualitative features of our hndings remain. can be constrained when 
~ 270°, because it predicts largest N{e). 5 can be weakly constrained when 5^“® = 


0° or 180°, it cannot be constrained when = 90°. The constraints on sin^26* 

’ MNS CHOOZ 

improves, but only slightly. 

Fig. 1 ^ shows the same constraints as in Fig. 0 but when the input data are calculated 
for sin^ 20*™® = 0.06. Now that the input sin^20 is chosen significantly below 

CHOOZ ^ CHOOZ ^ 


the GHOOZ bound eq. (|2.271) , the global picture of the constraint from the VLBL-I-LBL 
experiments with HIPA is more clearly seen. The constraints for the case with = 90° 
(top-right) and 270° (bottom-right) are now look more symmetric. In case of = 
90°, the same N{e) can be reproduced for the other angles by choosing appropriately 
small sin^20 , whereas for = 270°, the same N(e) can be obtained for larger 
sin^ The constraints for = 0° and 180° cases also look more or less symmetric. 

Locally, a small shift in can be compensated by a small shift in sin^ 20^^^^^^ to make 
N{e) similar at all experiments. This compensation does not work perfectly because the 
ratio V(e, 4GeV)/V(e, 6GeV) hits the minimum (~ 0.41) at around = 5°, and the 
maximum (~ 0.63) at around = 165°. These trends can be read off from Fig. 0- 
Unfortunately, the region allowed by < 4 extends over to all 5^j^g in both cases. From 
Figs. pTHl^ , we can tell that sin^ 20^^^^^^ can be constrained by these experiments, but 
the constraint can be improved significantly if can be constrained independently by 
other means. If indeed the LMA scenarios is chosen by the nature, it is important that we 
present the constraints on sin^ 20^^^^^^ and simultaneously. 
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^true _ no 

*^MNS ~ 

o 

o 

180° 

270° 

'\-min 

(A) 

1.0 

0 

0.7 

0 

(B) 

2.2 

0 

1.6 

0 

sin 20 (^jjqq 2 

(A) 

o.o50L;):i)(,^ 

o-048L);:()(,^ 

0.074L():()i^ 

o.o79L(j:;))^ 

(B) 

o.o59L;j:)();^ 

“oW^ 

0.068L(J:(Ji^ 

o.o80F():;j,)^ 

5m? 

ATM 

(xlO-3) (eV2) 

(A) 

o 4^5+0.056 

'J-^O-0.07 

q /lo+U.Uti 

'J-^^-0.09 

q qo+U.UY 

'J-'JO-0.085 

q A q+U.Ufci 

'J-^'J-0.09 

(B) 

q 7K;+u.ua 

■J- ' ^-0.05 

q /io+u.ut> 

'J-^^-0.04 

q qq+u ua 

'J-'jy-0.045 

Q A O+U.UDl^ 

^•^'^-0.045 

sin 20^.p^ 

(A) 

> 0.981 

> 0.975 

> 0.974 

> 0.979 

(B) 

> 0.992 

> 0.992 

> 0.993 

> 0.992 


Table 8: Constraints on the three-nentrino model parameters, siv? 
sin^26'^^j^, when SMA or VO scenario is assnmed in the £t. The inpnt data are generated 
for the same 4 points as those in Table (eq. ([4.26|) at sin^26'^™®^^ = 0.06 ). (A) Results 
with 500 kton-year each for NBB(T^peak = 6 GeV) and NBB(ii^peak = 4 GeV) at L = 2,100 
km. (B) In addition, 100 kton-year data from NBB((pjr) = 2 GeV) at L = 295 km are 
included in the £t. 


4.4.4 5m? and sin^ 29 

ATM ATM 

Finally, we study the capability of the VLBL experiment in measuring the atmospheric 
neutrino oscillation parameters, and sin^ 29^^^, accurately. Table shows the mean 

and one-sigma errors of 5m? and the one-sigma lower bound on sin^ when the data 
are calculated for 5m^ = 3.5 x 10“^eV^ and sin^26' = 1.0 in the LMA scenario with 

ATM ATM 

the hierarchy I. The input parameters are chosen for the LMA point of eq.( [4.26|) , but for 
sin^ 26**™® ^ = 0.06 and the four analysis. The £t has been performed by assuming 
the hierarchy I but allowing all the model parameters to vary freely within the LMA 
constraint eq. (p.32|) . The htting conditions are the same as in eq. (|4.29|) . The (A) rows 
gives the results with 500 kton-year each for NBB(i?peak = 6 GeV) and NBB(i?peak = 4 
GeV) at L = 2,100 km. The (B) rows gives the results when in addition, 100 kton-year 
data from NBB((p 7 r) = 2 GeV) at L = 295 km are included in the £t. The sensitivities to 
(5m^TM improved by using the data from the HIPA-to-SK experiments. 

The expected sensitivity for is about 3 ~ 4.5% and that for sin^ about 

2.4%. After including the HIPA-to-SK data the sensitivities improve to about 2 ~ 3% and 
1%, respectively. 

Table ^ shows the results of the £t to the same sets of data but when the SMA scenario 
is assumed in the £t; see eq. ([4.28| ). As we have seen in subsection 4.4.2 and in Fig. |^, 
the SMA (or LOW or VO) gives a reasonably good £t to the data that are calculated by 
assuming the LMA scenario. The minimum can become zero when = 90°, 270°, 
whereas is as low as 1 or 2 when <5*™® = 0° or 180°. On the other hand the fitted model 
parameters are signihcantly different from their input (‘true’) values. The magnitudes and 
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the significances of these differences between the inpnt parameters and their fitted valnes 
can have read off from Table & The shifts and errors of sin^ can be read off from 

Fig. |T^ Table ^ shows that the fitted can be more than 3% smaller than the inpnt 

valnes. It is therefore important that the SMA/LOW/VO scenarios are exclnded with good 
conhdence level, e.g. by KamLAND experiment |^, in order to constrain the parameters 
of the three-nentrino model by the LBL experiments. 


4.5 Results for L = 1,200 km 

In order to examine the sensitivity of the physics ontpnts of the VLBL experiments to the 
base-line length, we report the whole analysis for L = 1 , 200 km, which is approximately 
the distance between HIPA and Seonl. 

(a) NBB (Epeak=3GeV) at L=1200 km (b) NBB (Epeak=5GeV) at L=1200 km 
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Fignre 19: CP phase dependence of N{e) and N{^) at L = 1, 200 km for lOOkton-year with 
(a) NBB(i?pg 3 ^^ = 3 GeV) and (b)NBB(£'pg 3 ^]^ = 5 GeV). = 0° (solid-circle), 90° (solid- 
sqnare), 180° (open-circle), and 270° (open-sqnare). The inpnt parameters are = 

3.5 X lO-^eV^, sin^ 20^^^ = 1.0, 5ml^^ = 10 x lO-^eV^, sin^ 2 ^ 3 ^^ = 0.8, sin^ = 0.1, 

and p = 3 g/cm^. The predictions for the fonr types of the neutrino mass hierarchies (Fig. |I]) 
are depicted as I, II, III and IV. 


We show in Fig. [T^the expected signal event numbers, N{e) and V(/i), at the base-line 
length of L = 1, 200 km from HIPA with 500kton-year for (a) the NBB with F^peak = 3 
GeV and for (b) the NBB with Epeak = 5 GeV. The predictions are calculated for exactly 
the same three-neutrino model parameters and the matter density as in Fig. 0 and Fig. o 
; see eq. ( |4.16| ). The predictions for the neutrino mass hierarchy cases I to IV, see Fig. 
and Table |^, are shown by separate circles when the is allowed to vary freely. The 
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four representative phase values are shown by solid-circle (^^ns ~ 0°); solid-square (90°), 
open-circle (180°), and open-square (270°). The statistical errors of the N{e) and iV(/i) 
measurements at 500 kton-year are shown on the = 270° point for the hierarchy case 

I. 

Because we have learned from the analysis at L = 2,100 km the number of CC 
events, N{e), is most sensitive to the neutrino model parameters, we choose the peak 
energies (T^peak) of NBB by requiring large N{e) and suppressed iV(/i). A pair of 
NBB’s is then chosen such that the -dependence of the ratio of A^(e)’s is signihcant. 
The chosen peak energies, Ey = 3 GeV and 5 GeV at L = 1, 200 km have the same L/Ey 
with 77,^ ~ 5 GeV and 9 GeV at L = 2,100 km, respectively. As shown in the middle hgure 
in the left-hand-side in Fig. |T^, the magnitude of the transition probabilities for the 

mass hierarchy I is smaller than that for the hierarchy III in the interval 3 GeV ^Ey ^ 
6 GeV, whereas this ordering is reversed for 6 GeV ^Ey. Accordingly, in Fig. |^(b) for 
NBB(i7peak = 5 GeV) at L = 1,200 km, the predicted N{fi) in the hierarchy I is larger 
than that in the hierarchy III. It turns out that this reversing at the ordering of N{fi) 
is not very effective in distinguishing the neutrino hierarchy cases because the trend can 
easily be accounted for by shifting the atmospheric neutrino oscillation parameters, 
and sin^26*^^j^, slightly. The fact that the two chosen NBB’s cover the opposite sides of 
the Py^^y^ = 0 node leads to a signihcant improvement in the measurement of and 

sin^26' , see Table p|. 

When we compare Fig. ^ for L = 1, 200 km with Fig. for L = 2, 100 km, we notice 
that the reduction of N{e) for the hierarchy III at L = 1, 200 km is not as drastic as the 
reduction at L = 2,100 km, but that the expected statistical errors of the signal is smaller 
at L = 1, 200 km because of a factor of two larger N{e) for the same size of the detector. 

In Fig. 1^, we show the expected numbers of signal events, N{fi) and N{e), for the 
same set of NBB’s, (a) NBB (i7peak = 3GeV) and (b) NBB (i7peak = 5GeV), each with 
500 kton-year at L = 1, 200 km. The three-neutrino model parameters and the matter 
density used for calculating those numbers are the same as those used to generate Fig. 
for L = 2,100 km; see eq. (^4.17|) and eq. (|4.18|) . All the symbols are the same as those 
adopted in Fig. The dependences of N{e) on the input sin^ ‘^^caooz same as 

those found in Fig. |T^; N{e) decreases as input sin^ 26 *(.jjqq 2 is decreased from 0.1 to 0.08, 
0.06, 0.04, 0.02 and 0. The solid-circle, open-circle, solid-square and open-square points 
show the predictions of the LMA and SMA scenarios with neutrino mass hierarchy I, II, III 
and IV, respectively. For each hierarchy, the five larger grand circles give the predictions 
of the LMA scenario with = 15 x 10“® eV^, and the smaller circles are for the 

LMA with = 5 x 10“® eV^. The VO predictions of the neutrino mass hierarchies 

I and II (III and IV) cannot be distinguished, and they are given by the solid (open)- 
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(a) NBB( Ep33,=3GeV ) at L=1200 km 


(b) NBB( Epp^,=5GeV ) at L=1200 km 




Figure 20: The neutrino parameter dependences of the expected numbers of CC events 
and CC events, A^(e) and respectively, for the NBB with = 3 GeV (a) and 5 

GeV (b) 500 kton-year at L = 1, 200 km. All the symbols are the same as those in Fig. IB. 


diamonds. The difference in N{fi) is largest in the LMA scenarios with = 15 x 10“® 

eV^, for which the hierarchy III predicts about 40 % larger (20 % smaller) N{fi) than 
the predictions of the hierarchy I for the NBB with i^peak = 3 GeV (5 GeV). When we 
compare Fig. for L = 1,200 km with the corresponding Fig. for L = 2,100 km, 
we notice that the prediction for N{e) in the hierarchy III are signihcantly larger for 
L = 1, 200 km than those for L = 2 ,100 km. If it were only N{e) in Fig. ^ that effectively 
discriminates the neutrino mass hierarchy, we should expect signihcant reduction of the 
hierarchy discriminating capability of the VLBL experiments at L = 1, 200 km. 

In order to study these questions quantitatively, we repeat the fit for the following 
sets of experimental conditions ; 


(A’) 500 kton-year each for NBB(F^peak = 5 GeV) and NBB(F^peak = 3 GeV) 

at L = 1, 200 km 


(4.30a) 


(B’) In addition to (A’), 100 kton-year data from NBB((p^) = 2 GeV) 
at L = 295 km are included in the £t. 


(4.30b) 
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(a) Hierarchy I (true) vs III (fit) (b) LMA (true) vs SMA (fit) 



Figure 21: The same as Fig. |T5(b) and Fig. |T5|(b), but for L = 1, 200 km. 


In Fig. ^ we show the minimum as functions of the £t parameter sin^ ‘^^caooz 
data set (B’). The mean values of A^(/i, F^peak) and A^(e, i^peak) at F = 1, 200 km and those 
at F = 295 km are calculated for the experimental condition of eq. (|4.30b|) by assuming the 
LMA scenario and the hierarchy I, and by choosing 12 sets of the model parameters as 
in eq.(ra); sin2 2t = 0.02, 0.04, and 0.10 and = 0°, 90°, 180°, and 270°. The 


£t has then been performed by assuming the LMA scenario with the hierarchy III; see 


eq.( |4.27|) for details. The 12 lines in Fig. ^(a) correspond to sin 261^™®^^ = 0.02 (thin 
lines), 0.04 (medium-thick lines), 0.1 (thick lines) and <5^“® = 0° (solid lines), 90° (long- 
dashed lines), 180° (short-dashed lines), 270° (dot-dashed lines). We find that the neutrino 
mass hierarchy (between I and III) can be determined at 3a level, for all four values of 


i;™'atidwhensin^2»„„. 


= 0.1, and for = 180° and 270° when sin^ 26**™® = 0.04. 

’ MNS CHOOZ 


Amin 


is greater than about 7 for 5*™® = 0° and 90° when sin^ 20*^“® = 0.04. Even when 

^ MNS CHOOZ 


i 2 20true 

CHOOZ 


= 0.02, xLr, > 4 for = 180° and 270° 

7 /\-lfLLlL MNS 


These results at F = 1, 200 km 
are not much worse than the results at F = 2,100 km. We find that this is because the 
combination of the VLBL experiment and the F = 295 km LBL experiment is still effective 
in distinguishing the hierarchy cases, even though the VLBL experiment at F = 1, 200 km 


itself cannot distinguish the hierarchy cases when sin"^ 20"’“® ^0.04. 

CHOOZ 


In Fig. 1^ (b), the mean values of V(/i, Epeak) and V(e, Epeak) are calculated for the 
8 LMA points with hierarchy I; sin^ 20^^^^^ =0.06 and 0.10 for the four points in 
eq.( 4.26 ). The x^ fit has then been performed by assuming the SMA scenario within the 
hierarchy I; see eq.( 4.28 ). The results of the £t are given for sin^ ‘^^chooz ~ (thin lines), 
0.1 (thick lines) and = 0° (solid lines), 90° (long-dashed lines), 180° (short-dashed 
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Figure 22: The same as Fig. |T^, but for T = 1, 200 km. 


lines), 270° (dot-dashed lines). The qualitative features of the results shown in Fig. pT|(b) 
are very similar to those of Fig. [^(b) at L = 2,100 km. The VLBL experiment with 
HIPA do not seem to have the capability to distinguish the scenarios of the solar-neutrino 
oscillation models. This is mainly because the magnitude of oscillation phase A ^2 is smaller 
than 2.1 even for = 10 x 10“^eV^ and L/= 2,100km/4GeV. We should rely on the 

results of the solar-neutrino oscillation measurements and the forthcoming reactor neutrino 


oscillation experiments 

Fig. 


shows the regions in the sin^ '^•S- ^^mns pi^me which are allowed by the 

VLBL experiments at L = 2,100 km with 500 kton-year each for NBB (i7peak = 5GeV) 
and NBB (i7peak = 3GeV), and the LBL experiment at L = 295 km with 100 kton-year for 
NBB ((ptt) = 2GeV); eq. ([4.30b| ). All the model parameters used to calculate the expected 
numbers of events and the symbols of the hgures are the same as those adopted for Fig. p!8| . 
The four input points at V™® = 0°,90°, 180°, 270° and sin^ 26*^™® = 0.06 are shown by 
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5 =0° 

^MNS ^ 

90° 

180° 

270° 

drri'^ 

ATM 

(xlO-3) (eV2) 

(A’) 

q e:n+U.U65 

'J-JU o,o8 

'J-JU o.o8 

'J-JU o.o7 

(B’) 

q cfi+U.Ubt) 

'J-'-'U-O.OS 

'J-JU 0.075 

“q-ETTmUT 

'J-^U o.o75 

'J-JU o.o7 

sin"20^^^ 

(A’) 

> 0.9940 

> 0.9940 

> 0.9935 

> 0.9935 

(B’) 

> 0.9945 

> 0.9945 

> 0.9945 

> 0.9945 


Table 9: The mean and the la errors of 6 m] and the one-sigma level bound of sin^ 20*™® 
when the input data are calculated for the four LMA points (eq. (|4.26|) at sin^ ‘^^ckooz ~ 
0.06) and the hts are performed by assuming the LMA scenarios. (A’) Results with 500 
kton-year each for NBB(i?peak = 5 GeV) and NBB(i?peak = 3 GeV) at L = 1, 200 km. (B’) 
In addition, 100 kton-year data from NBB((p 7 r) = 2 GeV) at L = 295 km are included in 
the £t. 


the solid-circle in each plane, and the region with Xmin < and 9 are shown by the solid, 
dashed, and dotted boundaries, respectably. The fif has been performed by assuming 
the LMA scenario with the hierarchy I, but otherwise by allowing all the model parameters 
to vary freely; see eq. (|4.271) . 

The results shown in Fig. look very similar to those of Fig. ^ for the combination 
of the VLBL experiment at L = 2,100 km and the LBL experiment at L = 295 km; 
eq. ([4.25b|) . For = 90° and 270°, we cannot constrain the GP phase at all from these 
experiments. For = 0° and 180°, there appear a region of Xmin < 1 where is 
constrained, but the preferred region covers the whole range at Xmin < 4 level. We hnd 
that the area of the Xmin < ^ region is almost the same as that of the L = 2, 100 km 
experiment in Fig. for each case. We conclude that the capability of measuring 
sin^ 20 (^jjqq 2 in the LMA scenario is very similar between the VLBL experiments 

at L = 2,100 km and those at L = 1, 200 km. 

Finally, we show in Table ^ the expected accuracy of the measurements of 6771 ^^^^ and 
sin^20 . The expected number of events are calculated for = 3.5 x 10“^eV^ 

and sin^20^™^ = 1.0 in the LMA scenario with the hierarchy 1. The 4 LMA input points 
and the fit parameters are the same as those used for Table |^. The mean and the one- 
sigma errors for htted and the one-sigma lower bound for sin^ 29 are shown in 

the Table p. The results of the VLBL experiment only, eq. (|4.30aj) , are shown in the (A’) 
row, whereas those obtained by adding the LBL experiment at L = 295 km, eq. (|4.30b|) , 
are shown in the (B’) rows. By comparing the results of Table and Table we hnd 
that more precise measurements of and sin^ 20^.^.^ can be achieved by the VLBL 

experiments at L = 1, 200 km. As mentioned above, this is essentially because the two 
NBB’s that we choose for studying the potential of the L = 1, 200 km experiments cover 
the two sides of the = 0 node, whereas one of the two NBB’s chosen for studying 

the L = 2,100 km case sits on top of the node to minimize A(/i). We hnd that the 
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accuracy of and sin^ ‘^^caooz measurements at the L = 2 ,100 km experiments can 

be improved, if a NBB with higher peak energy is added. Because higher i^peak NBB gives 
higher r-backgrounds and larger N{fi), the optimal choice of the beam-type and beam 
energy should be determined according to the prime physics objectives and the details of 
the proposed detectors. 


5 Summary and discussions 


The KEK-JAERI joint project on HIPA (High Intensity Proton Accelerator) will start 
delivering 10^^ POT (protons on target) at 50 GeV per one operation year by the year 2007. 
This compares with 3.85 x 10^® POT at 12 GeV which are being provided by the KEK 
PS for the K2K neutrino oscillation experiment. The HIPA project thus provides us with an 
excellent opportunities of extending the K2K experiment, by shooting the neutrino beam 
from the HIPA to SK (the Super-Kamiokande) with L = 295 km |^. Rehnement of the 
measurements of the atmospheric neutrino oscillation parameters and sin^ ) 

and the hrst measurement of the —>■ Ve oscillation parameter (sin^ 26 '(^jjqq 2 ) by observing 
the electron appearance are the main targets of the proposed experiment m- 

In this report, we explore the possibility of extending the LBL experiments to a much 
longer base-line length (L > 1000 km) as a second stage of the neutrino oscillation exper¬ 
iments using the HIPA beam. In particular, we examine in detail the physics potential 
of VLBL experiments at L = 2,100 km, which is approximately the distance between the 
HIPA and Beijing, where strong interests in constructing a huge neutrino detector (BAND) 
have been expressed We also examine the case at L = 1,200 km carefully in order to 
study the sensitivity of physics outputs on the base-line lengths of the same order. The 
distance 1,200 km is approximately that between the HIPA and the central Korea, where 
strong interests in LBL neutrino oscillation experiments are expressed. 

We study the physics potential of such VLBL experiments within the three-neutrino 
model. The three-neutrino model gives a consistent picture of all the neutrino oscillation 
observations except the LSND experiment, and has six parameters those are observables 
by neutrino-oscillation phenomena: two mass-squared differences, 

(5m?n = ±(5m^ , three angles, sin^ 20 ^ , sin^ 20 . , and sin^ 29^ , and one phase, 6 

The main objectives of such VLBL experiments may be summarized by the following four 
questions : 


1. Gan we distinguish the neutrino mass hierarchy cases ? 

2. Gan we distinguish the solar neutrino oscillation scenarios and sin^20gQj^) ? 

3. Gan we measure the two unknown parameters of the model, sin^ ‘^^caooz ? 
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4. How much can we improve the measurements of and sin^ 

In order to quantify our answers to the above questions, we make the following simplihca- 
tions for the neutrino beams and the detector capability. 

(i) High-energy NBB’s with i^peak = 3 ~ 6 GeV for the VLBL experiments. 

(ii) Low-energy NBB with = 2 GeV to represent the LBL experiment to SK. 

(iii) 100 kton-level detector which is capable of detecting GG and z/g GG events almost 
perfectly, but is not necessarily capable of measuring the particle charges and hadron 
energies. 

(iv) Backgrounds from secondary beams and from leptonic r decays of the GG events, 
as well as those from ir^-e misidentihcation the neutral current processes are ac¬ 
counted for. 

(v) Gommon flux normalization errors of 3 % is assigned for the high-energy NBB’s, and 
an independent 3 % error for the low-energy NBB. 

(vi) Overall uncertainty in the matter density is assumed to be 3.3 % for each experiment. 

The detector capabilities of (iii) and (iv) are taken from those of the SK detector which 
may be common for water-Cerenkov detectors. 

We hnd the following observations. By assuming that the HIPA-to-SK experiment 
measures the atmospheric-neutrino oscillation parameters to be centered at = 3.5 x 

10 “^ eV^ and sin^ ^y assuming that its cumulative effects can be represented 

by 

100 kton-year at L = 295 km with NBB((p 7 r) = 2 GeV), 
we hnd that VLBL experiments of 

500 kton-year at L = 2,100 km with NBB(ii^peak = 6 GeV) 1 

500 kton-year at L = 2,100 km with NBB(ii^peak = 4 GeV) J 

can give the following answers : 

1. If the neutrino mass hierarchy I is realized in Nature, then the hierarchy HI can be 

rejected at 3a (la) level if > 0.04 (0.02). See Fig. (b). 

2. If the LMA scenario = 10 x lO”"^ eV^, sin^26'3Qj^ = 0.8) is realized in Nature, 

and if sin^ 26 *(.jjqq 2 >0.06, the SMA/VO scenarios can be rejected at one-sigma level 

when is around 0° or 180° but not at all when is around 90° or 270°. See 

Fig. ig (b). 


(5.1) 


(5.2) 
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3. If the LMA scenario = 10 x 10“® eV^, sin^26'gQj^ = 0.8) is realized in Nature, 

and if sin^ = 0.06, then 0.04 < sin^ 26 *(.jjqq 2 < 0-1 is obtained when 

is around 0° or 180°, 0.02 < sin^ 26 '(,jjqq 2 < 0-08 when (5j^j^g is around 90°, 0.045 < 
sin^ 26 *(.jjqq 2 < 0-12 when is around 270°. can be constrained to local values 
at one-sigma level when the true is around 0° or 180° but it is unconstrained 
when it is around 90° or 270°. See Fig. |^. 

4. If the LMA scenario = 10 x 10“^ eV^, sin^26'„_, = 0.8) is realized in Nature, 

i® measured with the 3% accuracy and sin^ ‘^^atm i'® level. See Table 

Summing up, a combination of LBL experiments at L = 295 km and VLBL experiments 
at L = 2,100 km with the NBB’s from HIPA can determine the neutrino mass hierarchy 
at 3a level if sin^ ‘^^caooz >0.04, cannot distinguish the LMA scenario from the SMA/VO 
scenarios well, can constrain sin^ ‘^^caooz ) ^^0 it can also constrain (5j^j^g at one-sigma level 
in some cases. The errors of the atmospheric-neutrino oscillation parameters are reduced 
to = ±0.1 X 10“^ eV^ and sin^ > 0.991. Very similar results are found for 

a combination of L = 295 km and L = 1, 200 km experiments. See Fig. Fig. and 
Table ||. 
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